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ABSTRACT 


Model test results were used to define the failure mechanism associated 

with the static penetration resistance of cohesionless and low-cohesion 

soils. Knowledge of this mechanism has permitted the development of a new 

analytical method for calculating the ultimate penetration resistance which 

explicitly accounts for penetrometer base apex angle and roughness, soil 

friction angle, and the ratio of penetration depth to base width. Curves 

relating the bearing capacity factors N and N to the soil friction angle 

c yq 

are presented for failure in general shear. 

Strength parameters and penetrometer interaction properties of a fine 
sand were determined and used as the basis for prediction of the penetration 
resistance encountered by wedge, cone, and flat-ended penetrometers of 
different surface roughness using the proposed analytical method. Because 
of the close agreement between predicted values and values measured in 
laboratory tests, it appears possible to deduce in-situ soil strength para- 
meters and their variation with depth from the results of static penetration 
tests . 

A procedure for determining the soil cohesion and friction angle from 
the results of static penetration tests is proposed. This procedure is 
illustrated by application to model test results, to penetration data 
presented by other investigators, and to penetration data obtained for the 
lunar surface by the Apollo 15 self-recording penetrometer and the Soviet 
Lunar Rover Lunokhod- 1. 
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PREFACE 


This report presents the results of one phase of the research under- 
taken as a part of National Aeronautics and Space Administration Grant 
NGR-05-003-406, "Lunar Soil Properties and Soil Mechanics," The in-depth 
study of the static penetration resistance of soils described herein was 
made as a background study in support of the Apollo Soil Mechanics Experi- 
ment (S-200) to aid in the analysis and interpretation of penetration data 
obtained on the lunar surface. 

The major part of the material presented in this report was developed 
by H. T. Durgunoglu for a dissertation in partial fulfillment of the 
requirements of the degree, Doctor of Philosophy in Geotechnical Engineering 
at the University of California, Berkeley. 
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CHAPTER ONE 


INTRODUCTION 

It has been recognized since early in the lunar exploration program 
that a knowledge of lunar soil properties is essential for satisfactory 
solution of a number of scientific and engineering problems. As examples 
are questions related to (1) formation and compaction of surface layers, 

(2) characterization of different materials, (3) slope stability and down- 
slope movement of soil and rock, (4) prediction of seismic velocities and 
analysis of seismic data, (5) estimation of thermal behavior, (6) character- 
ization of dielectric properties for use in radar and electrical studies, 

(7) gas diffusion through the lunar surface, (8) strength and deformation 
properties for use in engineering studies, and (9) definition of conditions 
for terrestrial simulation studies. 

To help in the development of this knowledge a soil mechanics investi- 
gation was included as a part of Apollo missions 11 and 12, and a formal 
Soil Mechanics Experiment was assigned to missions 14 through 17. As a 
part of this experiment it was desired that quantitative data be obtained to 
the extent possible within the constraints of lead time for apparatus 
development, spacecraft payload, and astronaut time and capability on the 
lunar surface. A self-recording penetration test device was selected for 
this purpose and included as a part of the Apollo 15 and Apollo 16 missions. 
A penetration test was selected because (1) it is simple, (2) it can be done 
by hand, (3) conditions can be explored to depths of several tens of centi- 
meters, and (4) there is much available terrestrial experience with this 
type of test. 

At the same time a number of limitations in the understanding of 
penetration resistance of soils in relation to soil strength and compressi- 
bility properties were recognized. Thus a basic investigation of pene- 
trometer-soil interaction was undertaken for the development of a suitable 
analytical technique for expression of the static penetration resistance of 
soils in terms of penetrometer base apex angle, base roughness, depth of 
penetration, and soil cohesion and friction angle. The specific objectives 
were: 
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(1) To develop an improved understanding of the soil failure mode 
during a static penetration test, 

(2) To determine the influence on penetration resistance of pene- 
trometer base apex angle, size, and roughness, relative depth of 
base, and angle of internal friction and compressibility of soil, 

(3) To develop an analytical technique which will allow prediction of 
static penetration resistance of soils in terms of physical 
properties of the penetrometer and basic soil properties or, 
alternatively, determination of in-situ soil properties from the 
results of penetration tests, 

(4) To explain previous test results, many of which cannot be 
adequately accounted for in terms of existing theories. 

This report presents the results of these studies. Initially, the 
state-of-the-art was critically reviewed (Chapter Two) . The failure 
mechanism associated with deep static penetration in soils was determined 
using model tests (Chapter Three). Knowledge of this mechanism permitted 
development of a new theory for determining the ultimate base resistance of 
a penetrometer in terms of factors that depend on soil friction angle, 
penetrometer base apex angle, depth to base width ratio, and penetrometer to 
soil friction (Chapter Four) . The strength and penetrometer interaction 
properties of a fine sand were then studied (Chapter Five) and the results 
were used as a basis for prediction of resistance to penetration by wedges, 
cones and flat-ended penetrometers of different surface roughness. Labora- 
tory penetration tests were conducted to provide a basis for comparison 
with predictions. The effects of each variable, (base apex angle, base 
roughness, relative depth, soil friction angle, and soil compressibility) 
were further investigated both experimentally and theoretically (Chapter Six) . 

Finally, a procedure for deduction of in-situ soil properties from the 
results of static cone penetration tests is outlined. This procedure is 
illustrated by application to model tests, to penetration data obtained from 
the literature, and to penetration resistance data for the lunar surface 
obtained during the Apollo 15 and Siviet Rover Lunokhod-1 missions (Chapter 
Seven) . 
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CHAPTER TWO 


BEARING CAPACITY AND 
PENETRATION RESISTANCE THEORIES 


THEORETICAL CONSIDERATIONS 
General 

Most theories for ultimate bearing capacity are based on the fundamental 
formula of Prandtl, which is valid for the case of a rigid-plastic, incompres- 
sible, weightless material whose shear strength characteristics are given by 
the Mohr-Coulomb failure criteria: 

“ c + atancj) (2.1) 

where T ^ - shear strength 
c = soil cohesion 
O = normal stress 

<j) « angle of internal friction of soil 

The solution of Prandtl (1921) for the bearing capacity under a strip 
load on a rigid-plastic, incompressible and weightless material is: 

q = cN c (2.2) 

where 

N c = cot<j> [ e Trtan( l > tan 2 (| + |) - 1] (2.3) 

Reissner (1924) considered the effect of a surcharge, q, and concluded that 
the bearing capacity was increased by an amount qN^ where: 

N - .*■* tan 2 (f + |) (2.4) 

q 4 2 


It can be seen that the bearing capacity factors N^ 
following equation: 


and N 

q 


are related by the 


N = (N - 1) cotcj) (2.5) 

c q 

A widely used bearing capacity equation which considers soil cohesion, 

friction, and surcharge was first presented by Terzaghi (1943) by combining 

Equations 2.3 and 2.4 and adding a friction term (1/2“Y BN ), giving: 

s y 
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q £ - cN + 1/2 *y BN + qN 
c 's Y q 


( 2 . 6 ) 


where 


q^ = ultimate bearing capacity or unit resistance under an 

infinitely long foundation of width, B 
c = soil cohesion 

Y g = unit weight of soil 

q = surcharge 

N , N , N = primary bearing capacity factors, f (40 . 

c Y 9. 


Inclusion of soil weight considerably complicates the mathematical 
solution. To handle this difficulty, many approximate methods have been 
proposed. 


Numerical Techniques 

The finite difference approximation based on the method of characteris- 
tics is widely used in the numerical analysis of differential equations of 
the type applicable to the bearing capacity problem. Lundgren and Mortensen 
(1953) used this method to obtain a solution of the strip footing bearing 
capacity problem. The same method was used for the axially symmetric 
circular footing problem by Cox et al (1961) , Cox (1962) , and Larkin (1968) . 

In recent years, the same technique has been applied to deep foundations. 
Strip loading conditions were analyzed by Graham (1968) and solutions were 
provided for axisymmetric circular foundations by Nowatzki (1971) and 
Nowatzki and Karafiath (1972) . 

Numerical methods have both advantages and disadvantages when compared 
to other methods for the calculation of bearing capacity. They may be 
summarized as follows: 

Advantages : 

(1) Soil weight can be considered easily and properly. 

(2) Dependence of soil friction angle on mean normal stress along a 
failure surface can be considered in the solution. 

Disadvantages : 

(1) The angle of internal friction and the unit weight of the soil are 

usually considered as independent variables, although, in fact, they 
may be interdependent. 
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(2) They are not conducive to development of simple graphs or 
formulations . 

(3) The Haar and von Karman hypotheses used in the analysis of axi- 
symmetric problems is not a good assumption (see Chapter Four) . 

(4) For deep foundations, the extensions of slip lines are arbitrary, 
e.g. Nowatzki (1971) has assumed that slip surface reverts to the 
penetrometer shaft as was assumed by Meyerhof (1951) (Fig. 2.1b); 
however, Nowatzki and Karafiath (1972) in the analysis of the same 
problem assumed that slip surface ends at the base level, which is 
equivalent to neglecting the shear strength of overburden, as was 
assumed by Terzaghi (1943) (Fig. 2.1a) . 

(5) Statically correct solutions found by the method of characteristics 
are kinematically inadmissible . 


Consequently kinematically more acceptable slip surfaces were introduced, 
e.g, Jumikis (1961), Balia (1962), Hu (1964), and Gorbunov-Possadov (1965). 

As yet, however, there is no completely correct theoretical solution for 
the friction (N y ) factor for the simple shallow footing problem. In view of 
this, Meyerhof (1955) and Brinch Hansen (1961) suggested the use of following 


empirical formulas for the calculation of the bearing capacity factor N v » 
Meyerhof 

Brinch Hansen N y = l„8(N q - 1) tancj) (2.8) 


N y = (N q - 1) tan (1.4(f)) 


Y 

(2.7) 


Slip Surface Assumptions for Shallow Foundations 

The bearing capacity equation (Equation 2.6) generally includes three 

bearing capacity factors (N , N and N ) . Commonly, the values of N and N 

c Y q c q 

are those corresponding to the slip surface valid for weightless soil 

<i.e.» q/(Y s B)=°°) . Further, the value of N y used is that corresponding to a 

slip surface valid for q/(Y B) « 0 (footing at the surface). 

s 

This situation leads to calculated ultimate bearing capacity values (in 

which Ny and are supposed to be determined independently) which are somewhat 

on the conservative side* This is due to the linear superposition of two 

terms, each of which is determined from a statically admissible rupture surface 

which is different from the surface which actually applies for values of 

q/(y B) other than zero or infinity* 
s 

If there is no surcharge other than the overburden, q/(y B) becomes equal 

s 

to D/B* In fact, when q/(y B) is greater than zero and less than infinity, 

3 
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(a) (b) 

Terzaghi Meyerhof 

FIG. 2.1 DIFFERENT FAILURE SURFACE ASSUMPTIONS 
FOR DEEP FOUNDATIONS 
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the slip surface is unique and intermediate to the two limiting slip surfaces 
shown at the top and bottom of Fig, 2.2. 

DeBeer (1967) has suggested the following correction formula for the 
bearing capacity calculated assuming two different slip surfaces: 

q f = V(l/2Y g BN Y + qN q ) (2,9) 

where and N q are bearing capacity factors calculated for two different 

slip surfaces and V is a correction factor that depends on <J> and D/B. The 

correction factor for <p = 30° is shown in Fig. 2.3a. Its maximum value is 

1.17 and occurs for q/(y B) or D/B equal to 0.4. Hansen and Christensen 

s 

(1969) calculated the same correction factor for a smooth strip footing, with 
the results shown in Fig* 2 0 3b* 

General Bearing Capacity Equation 

In Equation 2 0 6 the term (q N^) represents the increase in bearing 
capacity caused by overburden pressure* This term does not consider the 
effect of the shear strength of soil above the foundation level if the primary 
bearing capacity factors (which are valid for D/B = 0) are used* Therefore, 
for deep foundations new bearing capacity factors should be calculated, 
Meyerhof (1951) ® However, Brinch Hansen (1961) has suggested that this effect 
can be taken into account by depth factors (d , d , d ) . Also, because the 

c y q 

general bearing capacity equation is formulated for strip foundations, shape 

factors (£ , £ , £ ) must be employed when considering other foundation con- 
c T q 

figurations. With these modifying factors Equation 2.6 can be written in the 
following form: 


q. => cN K d + 1/2 y BN U + y DN % d (2.10) 

f c e e 's y y y s q^q q v ' 

where £ , E , £ = shape factors for the cohesion, friction* and 

c Y q 

surcharge terms. 

d , d . d = depth factors for the cohesion, friction, and 
c y q ' 

surcharge terms. 


Values of the primary bearing capacity factors, N^, N , and N 
Fig. 2.4. 


are given in 


Depth Factors 

Skempton (1951) proposed a depth factor for the cohesion term as 
follows : 
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FIG. 2.2 FAILURE PATTERNS UNDER A SHALLOW FOUNDATION 

{After De Beer, 1967) 
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A SHALLOW FOUNDATION {After De Beer, 1967) 



FIG. 2.3b CORRECTION FACTOR V ON LINEAR SUPERPOSITION 
ASSUMED BY BEARING CAPACITY FORMULA AS 
FUNCTION OF B/(B+D) FOUND FOR SMOOTH STRIP 
FOOTINGS (After Hansen and Christensen, 1969} 
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Primary bearing capacity factors 



FIG. 2.4 PRIMARY BEARING CAPACITY FACTORS 
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d c = 1.0 + 0.2 (D/B) (2.11) 

for purely cohesive soils (<j) = 0) and values of D/B up to 2.5. A constant 
depth factor d c = 1.5 was suggested for depths greater than D = 2.5B. 

Brinch Hansen (1961, 1966) suggested the approximate relationship.: 

d = 1.0 + 0.35 (D/B) (2.12) 

c 

as an average value of d £ for D/B smaller than 1.0 and values of <(> between 
0 and 40 degrees. Brinch Hansen (1961) has also proposed the following 
generalized and semiempirical equations for the depth factors: 


d 

e 


1.0 + 


0.35 


(B/D) + [0.6/ (1+7 tan **<(>) ] 


d Y 


1.0 


d 

q 




(2.13) 


These relations are shown in Fig. 2.5, The factors d £ and d^ approach 
the indicated limiting values with increasing relative depth, D/B for a 
given value of cj>. 

Meyerhof (1963) proposed the following equations for the depth factors 


d - 1.0 + 0.2 (D/B) tan (J + |) 


(2.14) 


for <t> = 0°: 


d - d =1.0 

q y 


for <j) > 10° and D/B <1.0: 


d q = d^ = 1.0 + 0.1 (D/B) tan (J + |) 


(2.15) 


The values of are slightly larger, whereas those of are smaller than 
the values proposed by Brinch Hansen. 

DeBeer (1967) suggested the following depth factors: 


d q - 1 + [tan 2 (45° - |)a ,,tan * - 1] e’™” <B/D) 


d y = 1.0 


(2.16) 
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FIG. 2.5 DEPTH FACTORS AS A FUNCTION 
OF RELATIVE DEPTH 
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Ultimate depth factors are estimated using the bearing capacity factors 
given by Meyerhof (1951) for shallow and deep foundations. A comparison with 
Brinch Hansen's ultimate values is given in Table 2.1, where it can be seen 
that Meyerhof's values are much higher. This is due to the nature of the 
failure mechanism assumed by Meyerhof (1951) for deep foundations, (see 
Fig. 2.1b), It should be noted that Biarez et al (1961) showed * from the 
results of model tests, that the slip surface reaches vertical tangency 
provided that the relative depth is sufficiently large. 


Table 2.1 ULTIMATE DEPTH FACTORS FOR COHESION TERM 


Angle of internal 
friction, ^-degrees 

Ultimate depth factor 
for cohesion term, (d c ) u i t< 

Meyerhof 

Brinch Hansen 

20 

3.3 

1.6 

25 

4.8 

1.7 

30 

6.3 

2.0 

35 

8.7 

2.6 

40 

14.4 

3.5 

45 

23.4 

5.7 

50 

36.2 

9.9 


Shape Factors 

Early suggestions by Terzaghi (1943) and Skempton (1951) concerning 
the shape factors may be summarized as follows: 

Circular Areas Square and Rectangular Areas 

5-1.3 5=1 + 0.2(B/L) 

c c 

5 y =0.6 5 y = 1 - 0.2(B/L) 

where B is the width and L is the length of contact area. 

Meyerhof (1951) presented a diagram for the determination of combined 
shape and depth factors that are functions of B/L, D/B, and <J>. Based on 
Meyerhof's diagram Brinch Hansen (1961) developed the following semi -empirical 
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equations for the individual shape factors that are independent of depth: 


5 = 1 + (0.2+tan 6 <f>) (B/L) 



% = 1 " f (0 • 2+tan 6 <f>) (B/L) 


(2.17) 


Meyerhof (1961b, 1963) also proposed the following empirical expressions for 
the individual shape factors: 


£ =1+0.2 tan 2 (45° + |) (B/L) 

c / 

£ = £ =1.0 for 4>=0° (2.18) 

q Y 

£ = £ = 1 + 0.1 tan 2 (45° + £) (B/L) for <j»10° 

q y l 


DeBeer (1967) presented the results of tests for determination of shape 
factor £^,. From these test results it was concluded that the shape factor 
£ has a value of 0.6 independent of the soil density. In a recent study, 
DeBeer (1970) proposed the following empirical formula for the shape factor 


V 


£ = 1 + (B/L) sin<}> 

q 


(2.19) 

Further evaluation of the proposed shape factors is made in Chapter Six. 


ADDITIONAL FACTORS INFLUENCING BEARING CAPACITY 

Dependence of Angle of Internal Friction on Mean Normal Stress 

According to Mohr-Coulomb yield criteria, soil strength may be char- 
acterized by Equation 2.1. However, Terzaghi (1925) and others have pointed 
out that the angle of internal friction (c}>) , for sand varies not only with 
density but also with the mean normal stress (cr^) , at a given density. Due to 
the fact that along a possible sliding surface underneath a foundation the 
values of normal stress are variable from point to point, the strength char- 
acteristics of a sand are not likely to conform to the simple linear function 
indicated by the Mohr-Coulomb yield criteria. 
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As listed under the advantages of numerical methods, it is more difficult 
to use a non-linear function than the linear Mohr-Coulomb equation. 

Yareshenko (1964) has suggested the following non-linear relationship be used 
instead of the Mohr-Coulomb criteria: 

T = (kcr ) 1/n (2.20) 

in 

where k and n are constants and functions of the angle of internal friction. 
The yield function given by Equation 2.20 has been employed in the solution 
of a plane footing problem by Berezantsev and Kovalev (1968). Later, Kingston 
and Spencer (1970) presented a solution procedure for a general non-linear 
function. 

An average value of mean normal stress along the slip surface is 
generally considered in analytical solutions* Meyerhof has shown that the 
average mean normal stress is approximately equal to one-tenth of the 
ultimate bearing capacity. The symbol ^ will be used to denote the mean 
value of mean normal stress along a slip surface. As shown in Fig. 2.6, the 
angle <f>' corresponding to the secant connecting the origin to the point 
CFg ^ was used (DeBeer, 1967) for the calculation of the ultimate bearing 
capacity. Thus the true curved strength envelope OMN Shown in Fig. 2.6 is 
replaced by the straight line OPQ intersecting the true curve at point P with 
an abscissa value of: 

CT g,M = q f /10 (2 * 21 > 

DeBeer (1967) gave the following empirical expression for the average mean 
normal stress along the slip surface: 

q f + 3q 

a g,M = ( 1 - sin< l > ) -^1 (2.22) 

The dependence of angle of internal friction on mean normal stress 
along a potential slip surface is of great importance in the interpretation 
of the results of model tests of shallow foundations. As the ultimate 
bearing capacity increases with the width of the footing for shallow founda- 
tions, so does the mean value of mean normal stress along the slip surface 
increase with the footing width, thus the secant angle <J>' decreases in value 
as may be seen by inspection of Fig. 2.6. As a result of this influence* 
the results of laboratory tests on very small shallow footings may lead to 
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Shear Stress, 


* 



Normal Stress , cr 


FIG. 2.6 DEPENDENCE OF <f> ON NORMAL STRESS 

(After De Beer, 1967) 
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an overestimation of the bearing capacity of much larger actual shallow 
footings. 

Progressive Rupture 

It may be inferred from the previous discussion of the influence of 
mean normal stress that it would be desirable to conduct loading tests on 
shallow footings with footings of actual full size. Such tests have been 
performed on a large scale by Muhs (1963) . The values of determined 
from these tests are shown in Fig. 2.7. When comparing these values to 
those determined in small scale tests, DeBeer (1967) points out that the 
large scale tests give larger values of N . than the small scale tests 
at low densities, while at high densities the reverse is true. 

This occurrence can. be explained by the phenomenon of progressive 
rupture. 

The progressive rupture phenomenon is due to the fact that during the 
gradual increase of load on a soil, the shear strength is not immediately 
mobilized at all points on the potential slip surface, but initially only at 
the points where the shearing stresses are largest. From these points, the 
rupture gradually extends to other points along the slip surface. This 
gradual progression causes modifications and variations of the soil properties 
along the slip surface. 

In loose soils, because of the compressional deformations which occur 
before the rupture load is reached, the density of the soil in the highly 
stressed zone is already increased before rupture, and the same is true for 
the shear strength, which increases with increased density. Therefore, at 
rupture, the shear strength corresponding to the initial density no longer 
governs, but there exists some variable shear strength along the rupture 
surface. 

In dense soils, the reverse occurs. In the highly stressed zones, the 
dense soil begins to dilate, causing a decrease in density and thus a 
decrease in shear strength to take place. Therefore, when the state of 
rupture along a slip surface in a dense soil is reached, the shear strength 
corresponding to the original density is not available along the entire 
rupture surface. 

DeBeer (1967) further states that the progressive rupture phenomenon 
is scale dependent for shallow footings because the relative settlement at 
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Bearing capacity factor. 


D,D 2 D 3 : Results of Muhs' (1963) tests. 

3 <f> determined from direct shear tests 

ABC: Mean curve of the experimental results. 

. 4> determined considering normal stress 
dependency 


AoBqC^ Mean curve of the experimental results. 
<f> determined from normal triaxial tests 
=19.6 to 78.4 kN/m 2 ) 



il i l i l i l 

20 25 30 35 40 45 50 

Angle of internal friction, <£- degrees 


FIG. 2.7 BEARING CAPACITY FACTOR N y DETERMINED 
FROM EXPERIMENTS (After De Beer, 1967) 
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rupture increases with the width of the footing. This dependency may be 
inferred from Fig. 2.8, which relates relative settlement at failure to 
footing size. 


Influence of Base Roughness 

The primary bearing capacity factors shown in Fig. 2.4 apply only to 
perfectly rough bases (6/ 4> = 1). Investigations by Meyerhof (1955) and 
others indicate that in cohesive soils the roughness has little 'influence 
on the bearing capacity. However, in cohesionless soils, the bearing 
capacity of a surface footing with a smooth base is significantly less 
than that for a footing with a rough base. To account for this difference, 
Meyerhof has suggested that the factor be multiplied by a roughness 
factor * (r ) , expressed as: 

r y = n r + 2 (1 “ n r 2) (2 ‘ 23) 


where n r is the degree of roughness, defined by the ratio of the tangents 
of the angle of base friction (6) , and the angle of internal friction (<j>) : 


n r “ tan6/tan<|> (2.24) 

Equation 2.23 applies only to plane, horizontal bases at the soil 
surface. 

Hansen and Christensen (1969) calculated N values for different 
values of base roughness and angle of internal friction. Their results are 
shown in Fig. 2.9, which indicates that the values of for perfectly smooth 
footings (6=0) are approximately one-half the values of for perfectly 
rough footings (6=4*) ♦ 

Recently, Graham and Stuart (1971) presented solutions showing the 
influence of base roughness and various assumptions as to the distribution of 
friction along the base upon N factors. A graphical summary of their 
solutions compared to other theoretical solutions is given in Fig. 2.10. 


Influence of Base Configuration 

The influences of various non-planar base configurations (e.g., wedges 
and cones) on bearing capacity factors have been obtained, with certain 
assumptions and under certain conditions, by Meyerhof (1961a). As an example, 
the bearing capacity factors for shallow strip loading on smooth and rough 
wedges in a soil with <j>=30° are shown in Fig. 2.11. 
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FIG. 2.8 RELATIVE SETTLEMENT AT FAILURE AS 
A FUNCTION OF FOUNDATION SIZE 

(After De Beer, 1967) 
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FJG. 2.9 BEARING CAPACITY FACTOR Hy FOR 

STRIP FOOTINGS AS A FUNCTION OF AND 8 
(After Hansen and Christensen, 1969) 
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30 32 34 36 38 40 42 44 

Angle of internal friction , <f> -degrees 


FIG. 2.10 COMPARISON OF THEORETICAL SOLUTIONS 
FOR N r FOR SURFACE FOOTINGS (Zero 
surcharge), (After Graham and Stuart, 1971 ) 
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Semiapex angle , a - degrees 


FIG. 2.H BEARING CAPACITY FACTORS FOR WEDGES 

(After Meyerhof, 1961a) 
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For perfectly rough wedges, there is little dependence of bearing 
capacity factors on the total apex angle (2a) for values of 2a greater 
than 90°. That is, for rough wedges and cones with obtuse apex angles, 
the bearing capacity factors are nearly equal to those for plane and hori- 
zontal contact areas. However, for perfectly smooth wedges, the values 

of N and N increase with increasing total apex angle, 
c q 

Values of N £ as a function of apex angle for perfectly smooth and 
rough cones at shallow and great depths in purely cohesive (<f> =0) soil 
as presented by Meyerhof (1961a) are shown in Fig. 2.12. 


Influence of Soil Compressibility 

Vesic (1967) has suggested that the relative compressibility of a sand 
mass may be expressed in terms of its rigidity index, I r » defined as: 


(1+V) (c+qtan<f>) 


(2.25) 


where E = elastic modulus 
c = soil cohesion 
<p = soil friction angle 
q = overburden pressure 
V = Poisson's ratio 

Bearing capacity factors calculated by Vesic (1967) using the assumption 
that the ultimate pressure on the soil cone under a foundation is equal to 
the ultimate pressure needed to expand a spherical cavity inside the same 
soil mass are given in Fig. 2.13. 

Vesic (1963) has also suggested that, for compressible soils, local (or 
punching) 1 shear failure, rather than general shear failure, occurs. Based 
on the shear pattern shown in Fig. 2.14a, the following expression for 
was developed: 

N = e 3,8 ^ taa<P tan 2 (45° + |) (2.26) 

q 2 

This equation is plotted and compared to the classic Reissner equation for 
11^ for general shear in Fig. 2.14b. It may be seen in Figs. 2.13 and 2.14 that 
for compressible soils (local shear conditions), the bearing capacity factor 
N is much lower than for incompressible soils (general shear conditions) . 
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Semi apex angle, a - degrees 


(a) SHALLOW DEPTHS 



Semiapex angle, a - degrees 


(b) GREAT DEPTHS 


FIG. 2.12 BEARING CAPACITY FACTOR N c 
FOR PURELY COHESIVE SOIL 

(After Meyerhof, I96lo) 
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6. 2.13 BEARING CAPACITY FACTOR N q 
FOR COMPRESSIBLE SOIL 

• (After Vesic, 1967) 
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FIG. 2.14a FAILURE MODE FOR LOCAL SHEAR 

(After Vesic, 1963) 



FIG. 2.14 b COMPARISON OF BEARING FACTOR N q 
FOR LOCAL AND GENERAL SHEAR 
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APPLICATIONS OF DIMENSIONAL ANALYSIS TO STATIC PENETRATION 


General 

Many problems may be analyzed by considering model tests that are 
assumed to give a true representation of prototype phenomena, but at a 
reduced scale. In order to apply classical dimensional analysis geometric, 
dynamic, and kinematic similarity requirements must be satisfied. 

The static penetration problem may include the following variables: 


= unit base resistance 
Y g = unit weight of soil 
B = base width of the penetrometer 
D = penetration depth 
<j> = soil friction angle 

c = soil cohesion 

6 = penetrometer to soil friction angle 

2a = penetrometer base apex angle (wedges or cones) 

Geometric similarity requires that 

'’'’’model ' ('’'’’prototype 
Dynamic similarity however, demands that: 


(2.27) 


’f'<V>model ’ ' '^s’’ prototype (2 - 28: 

Considerations of kinematic similarity may be ignored in the static 
penetration problem because the velocities involved are insignificant. 

Lundgren (1957) states that, in order to represent the static pene- 
tration test in dimensionally correct form, the following relationship 
should be used: 


D/B - f Cq f /(Y S B) ] (2.29) 

Hvorslev (1970) indicates that consideration of the dimensionless ratio 

q J (Y B) from the general bearing capacity equation is convenient in the 
I s 

investigation of data for different penetrometer sizes and soil conditions. 


Cohesive Soils (ij>=0) 

The bearing capacity equation for cohesive soils is normally expressed 
in the following form: 


q, = cN d + Y ON d 
H f c c 's 


q q 


(2.30) 
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For purely cohesive soils (<j>=0) , N * 1.0 and d = 1.0. Also, for D < 5B, and 

q . q — 

assuming that d £ is a linear function of relative depth (D/B) , Equation 
2.30 may be rewritten in dimensionless form: 

Y B 

q f /c = N c d c + (D/B) (2.31) 

For y g and c constant, but variable B, Equation 2.31 will yield families 
of curves in terms of the parameters q^/c and D/B. For shallow depths and 
small footing widths, the second term in Equation 2.31 is negligible and 
the equation may be written in the following form: 

q f /c = N c d c ‘ ' (2.32) 

which, represents a single curve for different sizes as shown in 
Fig. 2.15a. However, Equation 2.32 should not be used for high values of 
D and B or very small values of cohesion without first estimating the relative 
influence of the second term in Equation 2.31. , 

Cohesionless Soils (c a 0) 

The bearing capacity equation for cohesionless soils is normally 
expressed in the following form: 

% ‘ V, 1 , + ? 8*Vy - (2 ' 33 > 

As shown in the analysis of depth factors, d^ = 1.0 and for small depths 
(D<5B) is a linear function of relative depth (D/B) . Then Equation 

2.33 can. he written in dimensionless form: 

V<V> - y, £ + (2-W 

Therefore, a plot of the parameters q /y B vs. (D/B) yields a single curve 

1 S 

independent of the value of B, as shown in Fig. 2.15b. In recent years, 
various investigators have used these dimensionless parameters for graphical 
representation of the results of plate bearing tests on cohesionless soils. 
This equation does not, however, consider settlements that occur before the 
bearing capacity is developed. 

Cohesion-Friction (c-<|>) Soils 

The bearing capacity equation for c-<|> soils is normally expressed in 
the following form: 
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Bearing capacity ratio fii/yn) Bearing capacity ratio (^f/XB) 




0. CLAY b. SAND 

FIG. 2.15 BEARING CAPACITY RATIOS VERSUS RELATIVE 
DEPTH (After Hvorslev, 1970) 
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(2.35) 


q £ = cN d + y DN d + i y BN d 
H f c c 's q q 2 's Y Y 

For small depths, d £ and d^ are linear functions of relative depth (D/B) . The 

factor d^ equals 1.0 for all depths. Therefore, Equation 2.35 may be written 

in dimensionless form as follows: 

q f / = U/ (Y s B) 1M c d c + N q d q I + | N y < 2 • ■ ■ 

For c-<J> soils, the ratio q J (y B) is not independent of B, because c is a 

x s 

soil constant and the ratio c/(y B) decreases with increasing values of B. 

s 

That is, results of bearing capacity or penetration tests with penetrometers 
of various sizes will form a family of curves, and the values of q^/Cy B) for* 
a given value of D/B will decrease with increasing values of B. However, 
Equation 2.36 can be used to estimate c and <j> separately by having 
penetration data with two different sizes of cone. A proposed procedure for 
this calculation is outlined in Chapter Seven. 

STATIC PENETRATION TEST AND ITS APPLICATIONS 
Description 

Although static penetrometer equipment and procedures have not been 
standardized, some generally accepted practices have been developed. The 
Dutch Cone (Fig. 2.16a) is widely used. It has a base area of 10 sq cm and 
apex angle of 60°. The rate of penetration is from 15 to 20 mm/ sec. The 
normal sounding test (early version) has the following procedure: 

The tube is pushed together with rod and cone (see Fig. 2.16a) into 
soil for about 0.3 meters. Subsequently, only the inner rod is pushed 
downward while the tube is retained by the skin friction of the soil. The 
resistance to cone penetration is generally measured by a hydraulic gage. 

After penetration of the cone for about 125 mm, the tube is pushed down 
without pressure on the inner rod and cone. When the tubes are pressed 
further into the soil the cone moves down with the tube so that full 
frictional and cone resistance will be obtained. Tube and cone are 
advanced 75 mm after which the procedure is repeated thus giving measure- 
ments of total friction and cone resistance every 0.2 m. 

The sounding equipment was improved by Vermeiden (1948) , who developed 
the sleeve cone (see Fig. 2,16b) because of difficulties experienced with 
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BEGEM ANN FRICTION CONE PENETROMETER 

(1) Cone and friction sleeve retracted 

(2) Cone in extended position 

(3) Cone and friction sleeve both advanced 


the original cone associated with friction between tube and rod due to 
sand particles • In the procedure used by the Department of Municipal 
Works of Rotterdam, tube and rod are advanced simultaneously while the 
forces on tube and rod are measured separately. The measurements are 
generally performed automatically at the top of the rods by means of an 
electric pressure-gauge, 

Begemann (1965) reported the development of a friction cone (see 
Fig. 2.16c) enabling the measurement of local friction along a sleeve. 

The ratio of friction to cone resistance was found to be dependent on the 
type of soil, thus permitting an approximate determination of the soil 
profile. Later, electric strain gage penetrometers were developed, as 
described by De Ruiter (1971) . 


Determination of Soil Type 

Begemann (1965, 1969) has shown that there is a definite relationship 

between the ratio of unit frictional resistance (f ) to unit cone resistance 

s 


(q ) and the soil type as shown in Fig. 
c 

the following ratios: 

Soil Type 

soft rock or shells 
sand 
silt 
clay 


2.17, Schmertmann (1967) proposed 

V* c (%) 

0.0 - 0.5 
0.5 - 2.0 
2.0 - 5.0 
> 5.0 


Determination of Soil Compressibility 

There have been many attempts to relate compressibility of soils to 
the cone resistance. Bachelier and Parez (1965) gave the following relation- 
ship for oedometer modulus (E) , and cone resistance for silty and clayey soils 


E = 2.3 (q c /a) 

where a is a soil constant. 


(2.37) 


The values of a for different types of soil tested were in the range of 
0.33 to 1.0. 

Schmertmann (1970) proposed the following relationship between 
compressibility modulus of cohesionless soils and the cone resistance: 
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Cone resistance - kg/cm 



Percentage of soil particles 
smaller than I6u 





(2.38) 


E = 2.0 q 
s ^c 

where E = Young's modulus for sand in kg/cm 2 
s 

q c = Dutch cone bearing capacity in kg/cm 2 

Determination of Shear Strength of Cohesive Soils 

There have been many attempts to relate undrained shear strength of 
cohesive soils to cone penetration resistance. Gawith (1952) proposed the 
following relationships: 

c = q c /10 for soft clay (2.39) 

c = q^/14.8 for firm compact clay (2.40) 

where c = cohesion in kg/ cm 2 

q c = Dutch cone bearing capacity in kg/cm 2 

Later, Begemann (1965) gave the following relationship based on friction 
cone data: 

c = q c /14 (2.41) 

where c = cohesion in kg/cm 2 

q c = Dutch cone bearing capacity in kg/cm 2 


SUMMARY 

Primary bearing capacity factors based on classical Prandtl-Reissner 
solutions are widely used in Equation 2.6 to compute the ultimate bearing 
capacity of shallow foundations. The utilization of two different slip 
surfaces for the determination of the different bearing capacity factors 
(N c> N^, and N^) has been shown to underestimate the bearing capacity, in 
some cases by as much as 17 percent. 

The ultimate bearing capacity of deep foundations is generally 
determined using primary bearing capacity factors and ignoring the shear 
strength of the overburden as proposed by Terzaghi (1943). However, in 
general shear failure the shear strength of the overburden cannot be 
neglected. Brinch Hansen (1961) proposed modification of the primary 
bearing capacity factors with empirically determined depth factors for 
the determination of the unit bearing resistance of deep foundations. 

Many slip line geometries have been assumed for determining the 
bearing capacity of deep foundations, and the range of bearing capacity 
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values obtained from the different assumed failure mechanisms is large. 
Experimental justification for a proposed slip line geometry at failure 
for deep foundations has been provided by Biarez et al (1961) who showed 
that the slip surface reaches vertical tangency provided that the relative 
depth of the foundation is sufficiently large. The validity of the slip 
line geometry proposed by Biarez should be further investigated for the 
continuous deep penetration of penetrometers with different base configu- 
rations. Further, because of the large deformations involved, the 
fundamentals of the progressive development of failure surfaces during 
continuous penetration must be investigated. 

Meyerhof (1961a) investigated the influence of base configuration on 
the slip line geometry in the vicinity of the base. It was postulated 
that a soil cone (or wedge) forms under rough, blunt bases and that a plane 
shear zone forms adjacent to smooth, sharp bases. However, the validity 
of these assumptions needs experimental justification. Further, the 
influence of base configuration on slip line geometry in the vicinity of 
the base for intermediate base roughnesses (0<6/<J><l) should be studied. 

The extent of the plane shear zone, if any, should be formulated in terms 
of soil friction angle ($) and base roughness (<5/<|)) and the validity of 
any theoretical considerations should be justified by observations. 

The only analytical solution available for the determination of the 
bearing capacity of wedges (or cones) is given by Meyerhof (1961a) . 

This solution is valid only for very shallow and very deep foundations and 
perfectly rough or smooth bases. In this solution bearing capacity 
factors are given as a function of base apex angle and soil friction angle. 
The shortcomings of this solution can be summarized as follows: 

(1) Bearing capacity will be overestimated, even in general shear 
failure, due to the invalidity of the assumed failure mechanism 
for deep foundations (see Table 2.1). 

(2) For a given relative depth (D/B) of foundation, the bearing 
capacity can only be determined by linear interpolation between 
solutions for shallow (D/B<1) and deep foundations. Such inter- 
polation is later shown to be incorrect. 

(3) For a given intermediate base roughness (0<S/<jKl), the bearing 
capacity can only be determined by linear interpolation (later 
shown incorrect) between solutions corresponding to perfectly rough 
and perfectly smooth bases. 
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(4) No solutions are provided for determining N c when the soil 
friction angle is larger than 30°. 

Meyerhof (1955) and other investigators have indicated that for 
cohesive soils the roughness of a plane and horizontal contact area has 
little influence on bearing capacity. However, in cohesionless soils, 
the bearing capacity of a surface footing with a smooth base is only 
about half of that for a footing with a rough base^ 

Numerical solutions utilizing the Haar and von Karman hypothesis 
have been developed to determine the bearing capacity of circular shallow 
and deep foundations. However, these theoretical calculations are open 
to doubt, because the validity of Haar and von Karman hypothesis for soils 
is questionable. For this reason, it is concluded that theoretical calcu- 
lations of bearing capacities developed for circular foundations are still 
so uncertain that the best estimates are made by the application of 
empirical shape factors to the bearing capacity factors for strip founda- 
tions. Many empirical relationships for shape factors have been proposed, 
however, resulting in great range of values. Therefore, it is concluded 
that the validity of these relations over the range of interest should be 
investigated experimentally. 

Since general shear failure cannot reasonably be assumed for compres- 
sible soils, the influence of soil compressibility should be considered in 
determination of the bearing capacity (or unit penetration resistance) of 
foundations. There have been some empirical methods proposed which account 
for soil compressibility, but their general validity has not yet been 
established. 

For cohesionless soils, stress-strain relationships are stress, 
strain and density dependent. Because of progressive rupture, the proper 
selection of soil friction angles in the determination of the penetration 
resistance should be investigated. 

Many empirical relationships have been proposed to relate soil shear 
strength parameters, soil compressibility, and soil type to penetration 
resistance. However, no theory heretofore available has explicitly 
accounted for the influence of such important variables as base configura- 
tion, base roughness, and relative depth over the range of values of 
interest. The analytical procedures developed in this study which are 
presented in detail in the following chapters should provide a basis for 
overcoming some of these deficiencies » 
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CHAPTER THREE 


PAILURE MECHANISM ASSOCIATED 
WITH STATIC PENETRATION 


PREVIOUS STUDIES 

The failure mechanism associated with static penetration has been 
previously investigated both theoretically and experimentally. Terzaghi 
(1943) and Meyerhof (1951) both proposed that the slip line of the radial 
shear zone under a strip foundation (see Fig, 2.1) may be approximated by 
a logarithmic spiral (theoretically correct for weightless soils) . 

Meyerhof (1961a) reasoned theoretically that a plane shear zone exists 
adjacent to the penetrometer base, Biarez et al (1961) observed experi- 
mentally that a rigid wedge (or cone) develops under the base of rough 
flat-ended penetrometers for relative depths greater than one (D/B>1) . The 
fact that for great relative depths the slip line of the radial shear 
zone reaches vertical tangency was experimentally shown by Biarez et al 
(1961) and applied to pile foundations (flat-ended) by Hu (1965) . To 
date, no model studies have been reported which consider the effects of 
base configurations other than flat. Also, the effect of base roughness on 
the failure mechanism associated with wedge-shaped and conical penetrometers 
has not been studied. 

MODEL STUDIES 
General 

Model tests were performed to establish a rational basis for theoreti- 
cal developments and to improve the understanding of penetrometer-soil inter- 
action. The effects of penetrometer configuration, penetrometer to soil 
friction, soil relative density, and relative depth on the failure mechanism 
were investigated. 

Materials and Equipment 

The models were constructed in a special lucite box with inside 
dimensions of 15"xl5"xl5" (0.38m x 0.38m x 0.38m). The box consisted of 
three sections, each 5.0 inches (0.13m) wide, held together with tie bars 
and made watertight by a sealant tape. The box was so designed that the 
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tie bars can be removed and a section of box can be separated. A steel 
plate and electrical vibrator were used to densify samples. To prepare 
loose samples, the sand was deposited above a screen which was then slowly 
raised. In order to level the successive sand layers during placement a 
special leveling plate (screed) was used. During trimming, the sample box 
was placed on a tilted platform in order to safely trim the desired 
sections. A spatula and a sharpened trowel were used to trim the sample. 

The soil used in all the tests described herein was a fine air-dried 
sand, Monterey Sand No. 0 (engineering properties are given in Chapter 
'Five) . In order to facilitate observation of the failure mechanism, 
alternate horizontal layers of clean and colored sand were placed. The 
sand was colored using dye and carbon tetrachloride, a technique which 
minimizes the change in the characteristics of the sand. 

Model Preparation and Testing Procedure 

The models were prepared and tested as follows: 

(1) A filter paper was placed inside the box on top of the water 
exit, to serve as a drain and to prevent piping during later 
saturation and drainage of the samples. 

(2) The initial and subsequent alternate layers of clean and colored 
sand were placed and leveled carefully with the leveling plate. 
Changes in soil density due to placement of additional layers of 
sand were found to be insignificant. 

(3) The elevation of soil was measured and the average density was 
calculated. 

(4) Various penetrometers were pushed into the sand at a constant 
rate, 0.2 in. /min. (5 mm/min.) and the resistance was measured 
with a load cell and recorded using an electronic recorder print- 
out system. Penetration was measured by movement of a marker on 
a ruler. The test was stopped upon reaching the desired base 
elevation. 

(5) The penetrometer was disconnected from the load cell without dis- 
turbing the soil and left in the model. 

(6) The sand was then water-saturated slowly from the bottom under 
a very low hydraulic gradient. About four hours were required. 
The settlement due to saturation was negligible. 
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(7) The model was then allowed to drain overnight. The capillary 
tension remaining in the wet sand after drainage allowed the sand 
to be safely trimmed when the box was clamped to a wooden platform 
inclined approximately 35° to the horizontal, 

(8) The removable portion of the box was detached and the model trim- 
med along any desired section. Each section was photographed. 

Results and Observations 

More than 15 model tests were performed. The results of some of them 
are analyzed in the following paragraphs and summarized In Table 3,1. ' 

Photographs of the center sections of the models are given in Figs. 3.1 
through 3.8. Observations based on studies of the center sections may be 
summarized as follows: 

(1) For the tested range of relative depths (D/B) and relative densities, 
only general shear type failures were observed for tests using 
wedges. These test points are plotted on Fig. 3.9 which also shows 
the boundaries proposed by Vesic (1963) for different failure modes 
at various relative depths and relative densities. It may be seen 
that the limits of the general shear zone indicated by this research 
are wider than those proposed by Vesic for another sand. This 
difference is apparently due to the fact that the two sands do not 
have the same compressibility characteristics. 

(2) A plane shear zone exists adjacent to the penetrometer base and 
the topmost angle of this zone varies with the roughness of the 
base as may be seen in Figs. 3.4, 3.5, and 3.6. 

(3) Figs. 3.4, 3.5, and 3.6 also show a radial shear zone adjacent to 
the plane shear zone. It may be seen that the radial shear zone 
extends to the surface of the sand. Fig. 3.10 indicates that the 
slip surface of this zone may be closely approximated by a 
logarithmic spiral. 

(4) Figs. 3,4 shows that, for penetrometers whose base dimension 
exceeds the shaft dimension, a zone of loose soil is created 
along the shaft. 

(5) The approximate volume change of the soil during shear was calcu- 
lated by studying the cross-section photographs. It appears that 
the soil dilated during shear in each case. The dilation zones 
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Table 3.1 SUMMARY OF MODEL TESTS 
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FIG. 3.1a CENTER SECTION PHOTOGRAPH 
OF MODEL NO. i 


Unit base resistance, q f - psi 



FIG. 3.1 b PENETRATION RESISTANCE 

CURVE FOR MODEL NO. 1 
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OVERALL VIEW 



lb) CLOSE-UP VIEW 


FIG. 3.2 CENTER SECTION PHOTOGRAPHS 
OF MODEL NO. 2 
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(a) OVERALL VIEW 



(b) CLOSE-UP VIEW 


FIG. 3.3 CENTER SECTION PHOTOGRAPHS 
OF MODEL NO. 3 
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(a) OVERALL VIEW 



FIG. 3.5 CENTER SECTION PHOTOGRAPHS 
OF MODEL NO. 5 
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(a) OVERALL VIEW 



(b) CLOSE-UP VIEW 

FIG. 3.6 CENTER SECTION PHOTOGRAPHS 
OF MODEL NO. 6 
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CENTER SECTION PHOTOGRAPH 
OF MODEL NO. 8 
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FIG. 3.8a CENTER SECTION PHOTOGRAPH 
OF MODEL NO. 9 



FIG. 3.8 b CENTER SECTION PHOTOGRAPH 
OF MODEL NO. 10 
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Relative depth, 


Relative density , D r - % 



FIG. 3.9 TYPES OF FAILURE FOR FOUNDATIONS IN SAND 
AT DIFFERENT RELATIVE DEPTHS (%) 
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< t > , 430 

0.5 


FIG. 3.10 A COMPARISON OF OBSERVED AND PREDICTED FAILURE 
SURFACES FOR MODEL No. 5 
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were determined by measuring the distance increases between the 
dark-colored interfaces . 

(6) Very small rigid wedges were observed (Figs. 3.2 and 3.3) for the 
case of flat-ended penetrometers. This is believed (Abdul-Baki and 
Lewis, 1970) to be due to the fact that penetrometer to soil friction 
does not become fully developed for the flat-ended case. The low 
developed friction also causes the volume of the failed soil to 

be smaller than expected (see Table 3.2). 

(7) All of the photographs show that there are many shear surfaces 
developed during penetration. Each failure surface corresponds to 
a certain relative depth. The discontinuous, step-wise development 
of the shear surfaces is due to the fact that further penetration 
is required in order to develop the full shear resistance of the 
soil along another shear surface. 

(8) It may be seen from Figs. 3.7 and 3.8 that the soil zone undergoing 
shear is much smaller in the case of axisymmetric strain conditions 
than in plane strain. This is due to the fact that the deformations 
developed at any given section are much smaller than in plane 
strain. The three dimensional nature of axisymmetric strain condi- 
tions is also the probable cause of the fact that no distinct shear 
surfaces were observed for this case. 

(9) The affected volume of soil increases with an increase in the rough- 
ness of the penetrometer base (see Table 3.2). 

(10) The affected soil volume increases with a decrease in the base apex 
angle of the penetrometer (see Table 3.2 for Models 4 and 5). 

PROPOSED FAILURE MECHANISM FOR WEDGE LOADING 

A failure mechanism under strip loading by wedges which satisfies all of the 
observed features of the model tests as well as the knowledge obtained from pre- 
vious studies for the range of relative depths considered is shown graphically 
in Figs. 3.11a and 3.11b. A plane shear zone exists adjacent to the base of 
the penetrometer (see Fig. 3.11b). A logarithmic spiral approximates the slip 
surface of the radial shear zone and either intersects the ground surface at 
point E (see Fig. 3.11a) or becomes vertically tangent to line EF (see Fig. 
3.11b), A similar failure mechanism was proposed for flat-ended piles by Hu 
(1965) and by Abdul-Baki and Lewis (1970). 
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Table 3.2 COMPARISON OF OBSERVED AND PREDICTED GEOMETRIC CONFIGURATIONS 
OF FAILURE SURFACES 
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Referring to Fig* 3.11b for the more general deep penetration case, it 

has been shown by Abdul-Baki and Lewis (1970) that the point 0 is the center 

of the logarithmic spiral. The zone bounded by OCEFG is one in which the 

state of plastic equilibrium has been reached at every point. The earth 

pressures developed along vertical faces EF and 0G are considered. It is 

assumed, as proposed by Hu (1970) that no shear stress develops along the 

vertical face EF. This implies that for penetration depths greater than D 0 , 

P 

all the soil deformation is accommodated below plane HE. The effects of the 
angle of internal friqtion and the penetrometer to soil friction on the 
magnitude of the topmost angle of the plane shear zone (y) are considered in 
the following section. Further, for blunt, rough (6»<j>) bases, model tests 
show that a rigid soil wedge (or cone) will be developed in front of the tip 
having a base angle of 45°+<}>/2. 

A comparison of predicted and observed geometric features of the failure 
surfaces for the model tests is presented in Table 3.2. 

THEORETICAL CONSIDERATIONS 

Determination of the Topmost Angle (y) of the 
Plane Shear Zone 

The geometric configuration of the plane shear zone adjacent to the 

wedge (Fig. 3.11b) is determined by the known wedge semi -apex angle (a), the 

topmost angle (y) , and the included angle ACO which is equal to (90* - <j>) . 

As the roughness (6/<J)) of the wedge increases, the angle y at point 0 

decreases and vanishes for a perfectly rough (6=<j>) wedge. Means for calcu- 
lating this angle are given in the following paragraphs. 

Cohesionless soils (c=0) 

As shown in Fig. 3.11b, the logarithmic spiral starts from point C. 
According to plasticity theory OC should be a plane along which full mobiliza- 
tion of shear strength of soil takes place. In other words: 

x b = cr b tan<J> (3#1) 

in which x b = shear stress on plane OC, a = normal stress on plane OC, and 
<J> = angle of internal friction of soil. It should also be noted that the 
stresses on plane OA should satisfy the following relationship: 

T a = a a tan(S (3.2) 
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in which T = shear stress on plane OA, a = normal stress on plane OA, 

ci a 

and 6 = penetrometer to soil friction angle. 


It can be seen from Fig. 3.12a that the following relationships can be 
written: 


and 


x 

a 


, lf~ g 3f 
K 2 


cos ( 2 y- 4 >) 


a 

a 



+0 

2 ~ 




3f 


■) sin( 2 y-<j>) 


(3.3) 


where •=* major principal stress at failure, “ m i nor principal stress 
at failure and y ■ the topmost angle of the plane shear zone. By substituting 
Equation 3.3 into Equation 3.2, the following relationship may be obtained. 


tan 6 


(°lf _a 3 f> cos ( 2 Y-4>) 

(aif+a 3f ) + ( a lf - a 3f ) sin( 2 y~<j)) 


(3.4) 


or 


tan 6 


[(ai/a 3 )f-l] cos(2y-<J>) 

[ (^ 1 / 03 ) f+i] + [( a i/a 3 )f- 1 ] sin( 2 y-cj)) 


By introducing K = Equation 3.5 becomes: 


(3.5) 


tan<$ = 


cos( 2 y-(j>) 

[ (K+l) / (K-l) ] + sin(2y-<()) 


(3.6) 


From Fig. 3.12a, the ratio of major principal stress to minor principal 
stress at failure can be written as follows: 


K - (a 1 /a 3 ) f 


1 + sin(}> 
1 - sincj) 


By substituting Equation 3.7 into Equation 3.6: 


tan 6 = 


cos( 2 y-(j)) 


(l/sin<j>) + sin( 2 y-<j>) 


(3.7) 


(3.8) 


or 


tan 6 [l+sin<{> sin( 2 y-<j>)] - sin<j) cos( 2 y-(j>) « 0 (3.9) 
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FIG. 3.12 DETERMINATION OF THE TOPMOST 
ANGLE (r ) OF THE PLANE SHEAR 
ZONE 
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Equation 3.9 indicates that for cohesionless soils, there is a unique value 
of angle y, for given values of 6 and <|>. 

For perfectly smooth wedge (6=0) in cohesionless soils. Equation 3.9 
becomes: 

sincf) * cos(2y-(j)) = 0 (3.10) 

For #0, cos(2y-<j)) =0 or 2y-<}> = (2n+l) j and consequently (for n=l); 

y = 45° + | (3.11) 

This may be seen directly from Fig. 3.12a. For 6=0, the line OL coincides 
with the o axis and 2y becomes equal to 90°+(f>. 

For a perfectly rough wedge (6=<J>) , Equation 3.9 becomes: 


tan<j> [l+sin<J> sin(2y-<|>)] - sin<(> cos(2y-r<|>) = 0 (3.12a) 

or for (j>^0: 

1 + sin<j> sin(2y-<j)) - cos<f> cos(2y-(J>) ■ 0 (3.12b) 


or 

1 - cos(2y) = 0 


and consequently; 


y = 0 

This also can be seen from Fig. 3.12a. For 6=<f>, lines OL and OS coincide 
and the angle y vanishes. 


Cohesive soils (6=0) 

For cohesive soils (<J>=0) , Equation 3.8 does not apply. In this case, 
the angle y can be calculated as follows: The roughness of a wedge for 

4>=0 soils is defined as: 


f 


e 


c /c 
a 


(3.13) 


where 0<f <1 is 
— c— 

wedge and soil) , 


the roughness factor, c 
and c is the cohesion. 

sin (f - 2T) = ^ - -J 6 


is the adhesion (cohesion between 
From Fig. 3.12b: 

(3.14) 
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or 

Y = j ~ \ sin -1 (f c ) (3.15) 

For perfectly smooth wedges, = 0 and y becomes equal to 45°, and for 
perfectly rough wedges, f^ = 1 and y becomes equal to zero. 

Cohesion-friction (c— <j>) soils 

Similar to cohesionless soils, the shear stresses on plane OC are given 
by (see Fig. 3.13): 

= c + tan<j) (3.16) 

and the shear stresses on plane OA are given by: 

T = c + O tan <5 (3.17) 

a a a 

Also, from Fig. 3.13, T and a are given by Equation 3.3. By substituting 

cL a. 

Equation 3.3 into Equation 3.17: 

(tfl f -a 3 f) cos(2y-4>) - 2c & + [(a lf +a 3f )+(0 lf -a 3f )sin(2y-(j))]tan6 

(3.18) 


The major principal stress at failure is given by: 

_ 1 + sin(j) 2c cos(|) 

If 1 - sintj) 1 - sin<f> 

By substituting into Equation 3.18: 


(3.19) 


F (c a , c, 6 , <p, y, a 3f ) = 0 (3.20) 

which implies that for c-<p soils, the angle y is a function of can 

be seen, however, that for a perfectly rough wedge Y=0, and for a perfectly 
smooth wedge y - 45° + y . 

It can be shown that the dependency of the angle y on 0 vanishes under 
the following assumption: 

c /c = tan6/tancj> (3.21) 

a 

and that functional F_ becomes: 

F (c , c, 6, <j), Y) = 0 (3.22) 

a 
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Vc * ton 8 /ton <t> 


FIG. 3 13 DETERMINATION OF THE TOPMOST 
ANGLE (r) OF THE PLANE SHEAR 
ZONE FOR COHESION - FRICTION 
. (C- j SOILS 
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Equation 3.21 implies that lines TS and T'L intersect at point 0* on the 
axis (see Fig. 3.13). It can thus be seen that functional F may be 
calculated. 

Geometrically it can be seen that 00‘ - c cot<J> = c a cot6 or that 

c /c = tan<5/tan<j>, which is a reasonable assumption for most c-<j> soils. With 
a 

this assumption. Fig. 3.13 yields the following: 


MB » KB/cos<j> 

and AA^ =* MA cos(2y-<j)) 


(3.23) 


Since MA^MB^radius of Mohr's circle: 


AA 1 - KB (3.24) 

From Fig. 3.13: 

O'K = KB cotcj) 

KM = KB tan(j) (3 . 25) 

M^l “ AA^tan(2y~^) 

O f A^ - AA^ cot 6 

By substituting Equation 3.25 into O’A^O'K+KMH- MA^, the following 
relationship is obtained 

AA^ cotS = AA^ tan(2y-<{)) + KB(cotc}> 4* tan<j)) (3.26) 


By substituting 


Equation 3.24 into Equation 3.26: 

cot6 = tan(2y-(j)) H (cotcj) + tan(f>) 

cos(zy-cp; 


or, 


tan 6 fl + sin(f) sin(2y-(j)) ] - sincj) cos(2y-c})) = 0 (3.27) 


It can be seen that Equation 3.27 is exactly the same as Equation 3.9. 

In other words, the angle y satisfies the same relationship in both 
cohesionless and cohesion-friction soils provided that Equation 3.21 is 
satisfied. 

Values of the angle y for different values of roughness and angle of 
internal friction may be calculated from Equation 3.9 using Subroutine AN G 
(See Appendix A) . An iterative procedure was used to calculate values of 
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Y to an accuracy of 0.1° . These values are tabulated in Table 3.3. Fig. 3.14 
shows the variation of angle y with roughness for different values of angle 
of internal friction. 

The following conclusions can be drawn from the analysis of Table 3.3 
and Fig . 3.14: 

(1) The effect of angle of internal friction on the angle y is much 
more pronounced for low values of roughness than for high values. 

(2) For high values of roughness, there is a considerable increase in y 
values for a small decrease in 6/t)> values. 

(3) The theoretical curve for y versus <S/<j> is well above the linear 
variation curve, indicating that linear interpolation for y values 
between 0° and 45°+<f>/2 for different roughness values is not 
permissible. 

Determination of Critical Relative Depth 

The vertical tangency point of the shear surface coincides with the 
ground surface (see Fig. 3.11b) for a certain relative depth depending on 
base apex angle, base roughness, and angle of internal friction of the soil. 
This relative depth is defined as the critical relative depth, (D/B) cr . If 
the relative depth of the penetrometer base is greater than the critical 
relative depth, the angle 3 (see Fig. 3.11b) will be equal to the angle of 
internal friction (<{>); otherwise it will be smaller t:han the friction 
angle and must be calculated by iterative procedures. The iterative technique 
used is described in Chapter Four. A summary of calculated critical 
relative depth values is given in Table 3.4. The significance of these 
values can be summarized as follows: 

(1) If the relative depth is greater than the critical relative depth: 

(a) For further penetration, there will be no change in the 
bearing capacity cohesion factor (N £ ) values (see Chapter 
Four) . 

(b) For further penetration, the increase in friction- 
surcharge bearing capacity factor (N ) will be propor- 
tional to the increase in depth. In other words, N 
values for larger depths can be calculated by linear 
extrapolation (see Chapter Four) . 

(2) They are used as the basis for the calculation of depth factors (see 
Chapter Four) . 
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FIG. 3.14 EFFECT OF ROUGHNESS ON ANGLE 
FOR GIVEN VALUES OF ANGLE OF 
INTERNAL FRICTION 
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Table 3.4 VALUES OF CRITICAL RELATIVE DEPTH, (D/B) 


a=15° 

/<? 

0.0 

25 

1.37 

30 

2.06 

35 

3.10 

40 

4.75 

45 

7.55 

50 

12.72 


0.0 0.5 1.0 


a=45‘ 


0.0 0.5 1.0 


1.38 
2.27 
4.01 

3.06 I 4.95 I 7.38 
14.48 

50 i 8.99 117.34 131.34 


\ 6 / <p 

o 

• 

o 

0.5 

1.0 

25 

0.64 

0.89 

1.31 

30 

1.02 

1.48 

2.27 

35 

1.64 

2.49 

4.01 

40 

2.70 

4.36 

7.38 

45 

4.66 

| 

8.14 

14.48 

50 

8.68 

16.75 

31.34 























Determination of Lateral Extent of the Radial Shear Zone 


The lateral extent of the shear surface is a maximum for relative 
depths equal to or greater than the critical relative depth. The lateral 
distance from the pole of the logarithmic spiral to the shear surface at 
the point of vertical tangency is denoted by r^ (see Fig. 3.11b) . Table 3.5 
summarizes the maximum lateral extent of the slip surface, in terms of the 
dimensionless variable r^/B, for various values of base semiapex angle (a), 
base roughness (6 / <f>) , and soil friction angle (<f>) . 

The importance of considering the lateral extent of the failure surface 
has recently discussed by Nowatzki and Karafiath (1972) . The soil volume 
affected in the two dimensional (plane strain) problem will be proportional 
to r c and the volume affected in the three dimensional (axisymmetric) problem 
is proportional to r^ 2 . Volumetric considerations are particularly important 
in compressible soils because to develop full frictional resistance along a 
larger failure surface larger soil volume must be compressed. Consequently, 
indices obtained with penetrometers whose configurations result in large 
affected soil volume are likely to be more representative of soil compressi- 
bility than of soil shear strength properties . 

SUMMARY 

Previous studies of the failure mechanism associated with static 
penetration have been reviewed. Model tests have been performed and the 
results used to further define the mechanism of failure. The effects of 
penetrometer configuration, penetrometer to soil friction, soil friction angle, 
and relative depth on the failure mechanism have been investigated. A 
proposed failure mechanism for wedge penetration based oh the interpretation 
of features observed in the model tests and the results of previous studies 
is presented. The geometric features of the proposed failure mechanism have 
been formulated in terms of base semiapex angle, base roughness, soil friction 
angle, and relative depth of penetrometer 'base. 
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Table 3.5 SUMMARY OF r /B VALUES FOR D/B>.(D/B) 

G C2T 
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CHAPTER FOUR 


THEORETICAL CALCULATION OF STATIC 
PENETRATION RESISTANCE 


INTRODUCTION 

Both the results of model tests (Chapter Three) and additional experi- 
mental results (Chapter Six) show that for wedge (or cone) penetrometers, 
the ultimate penetration resistance depends strongly on the base semiapex 
angle (a), base roughness (<5/4>), and relative depth (D/B) . No theory avail- 
able heretofore has considered these three parameters explicitly over the 
ranges of interest. Knowledge of the failure mechanism has now permitted 
development of new relationships for the ultimate base resistance which 
account explicitly for these parameters. The development of the these 
relationships is described in this chapter, and curves showing bearing 
capacity factors versus angle of internal friction for several values of the 
parameters a, 6/<J>, and D/B are presented. The importance of various para- 
meters is discussed. 

GENERAL CONSIDERATIONS 

The general bearing capacity equation for a plane, horizontal strip 
foundation is: 


q f " cN c + \ T s BN y + qN q (4,1) 

where q^ = ultimate bearing capacity, q = surcharge, c = cohesion, y g = 

effective unit weight of soil, B = width of loaded area, and N , N , and 

c y 

N are bearing capacity factors for cohesion, friction and surcharge 

respectively. When the general bearing capacity equation (Equation 4.1) 

is applied to square, circular or any other contact areas of limited extent, 

the primary bearing capacity factors must be modified by shape factors 

(€.>€*€) which are usually determined empirically. With these modifying 
c y q 

factors, Equation 4.1 is expanded to: 

If ■ cN c«c + I Wy + «Vq (4 ' 2) 
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It is customary (Meyerhof, 1951) to calculate II and N for one assumed 

c q 

failure surface and N for another. As discussed in Chapter Two, the use 
of factors determined in this manner may lead to a significant underestimat- 
ion of the ultimate bearing capacity. As an alternative to the customary 
procedure, the factors N and N may be combined (N ) and a proper single 

H I 

failure surface may be considered for the calculation of and (Meyerhof, 

1951 and Hu, 1965). Consequently, Equation 4.2 may be rewritten as: 

q, = cN C + y BN £ (4.3) 

H f c^c Vs yq yq v 7 


where N is the beating capacity factor for the friction-surcharge term 
and is the corresponding shape factor . As previously mentioned, for 
wedge or cone penetrometers, both and are functions of the following 
variables : 


(1) Soil friction angle (<()) 

(2) Base semiapex angle (a) 

(3) Base roughness (<5/<j>) 

(4) Relative depth of penetrometer base (D/B) 

Thus : 


N c , N = f v f 2 (4>, a, 6/4), D/B) (4.4) 

DETERMINATION OF BEARING CAPACITY FACTORS 
General 

Fig. 4.1 shows the free body diagram for the calculation of the bearing 
capacity factor N c . The following expression for N^ may be derived (see 
Appendix A for details) : 

[l+sin4> sin(2y-4>)] _20 tan4> cos(2y-4>) tani|) n 20_tan4> 
c sinip coscp coscp 

+ [sin(2£+4>)-sin4)] [l+sin4) sin(2y~4>)] 20 tan4) _ 1 

cos 2 4>[cos4>-tan4>[sin(2£+4>)-sln4>] ] tan(p 

tan4> tan4> cos(2y~4>) [sin(2g+4))-sin4>] 20 tan4> ,, 

cos4>[cos4>-tan4>[sin(2£+4>)-sin4>] ] 6 ° 
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FIG. 4.1 FREE BODY DIAGRAM FOR DETERMINATION 
OF BEARING CAPACITY FACTOR N c 
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Because it is assumed (Hu, 1970) that there is no shear stress developed 
on the vertical plane CD, the angle £ vanishes, and Equation 4.5 simplifies 
to: 

vr _ ltsintj) sin(2Y-(j>) 2 0 tantj) 1 cos (2Y-4>)tan4) 2 0 tantj) 

c sintj) cos4> tan4> costj) e ° (4.6) 

where N c = bearing capacity factor for cohesion term, 

<j) = soil friction angle, 

y = the topmost angle of the plane shear zone, 
ip = 90°-rtx (a « semiapex angle) , 

9 o * 180°-(4H-y)+3 

Equation 4.6 can be written in functional form as: 

N c = F C4* > Y, 6) • (4.7) 

or because ^=90°-a, Y~^.( » <5 / ) and B=f^(<j>, 6/<j>, a, D/B), Equation 4.7 can 
be restated as: 

N c = f 1 (<j>, a, 6/4>, D/B) (4.8) 

Fig. 4.2 shows the free body diagram for the calculation of bearing 
capacity factor .N . From static equilibrium of body OCEFG: 


EM =0 
o 


(4.9a) 


or 


P b L b + P S 1 L ,1 ■ P s2 L ,2 + P l L pl + P 2 L p2 + W l L wl 


(4.9b) 


The following expression for N may be derived (see Appendix A for details) 


N = cosQ|)-6) [l+sin4) sin(2Y-4))] ) cos 2 (Y-<j)) T , 

Yq cos6 costj) eos(Y~<j>) \ 4cos 2 ij> cos 2 4> X 0 


3cos(y-4>) cos 2 8 2 0 tantj), 2 IA „ costj) costj) 
4co# costj) e ° (m ~ 3 m } " K o cos (Y-tj)) 


• (m-m ' ) • (m+2m ' ) + K 


COS 4) costj) 

cos(y~4 ) ) 



tantj) 

4 


(4.9) 
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FIG. 4.2 FREE BODY DIAGRAM FOR DETERMINATION 
OF BEARING CAPACITY FACTOR N yq 
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where 


yq = bearing capacity factor for friction-surcharge, 

$ = angle of internal friction of soil, 

ip = 90°-a (a=semiapex angle) , 

5 = base to soil friction angle, 

y = the topmost angle of the plane shear zone, 

e Q = 180° -(^+Y) +3, 

K = lateral earth pressure coefficient (see the discussion later 
in this chapter for proper selection) , 
m = relative depth (D/B) , 

»' = Dg/B, 

Dg = the vertical distance of point E on the failure surface 

above base level {a function of 8 (see Fig. 4.3)], 

t _ A sin3 cos(Y~ ( i ) ) 9 tan<p 
m 2 cosip cos<j> 6 ° 



(4.10) 


(4.11) 


Determination of the Angle 3 

In order to calculate the bearing capacity factors and N from 
Equations 4.6 and 4.9, the value of the angle -3 (see Fig. 4.1) must be known. 
For relative depths equal or greater than the critical relative depth, the 
angle 3 is equal to the angle of internal friction (<j>) of the soil. For 
relative depths less than the critical relative depth, the failure surface 
will intersect the ground surface before reaching vertical tangency. In 
this case, 3 will be smaller than cf> and must be calculated by iterative 
procedures. The procedure used for the calcualtion of 3 is described in the 
following paragraphs. 
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m 4.3 ITERATIVE PROCEDURE USED FOR THE 
CALCULATION OF CORRECT VALUE OF 
ANGLE 0 (j8 c ) 
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From Fig. 4.3: 


AB = B/ (2cos ip) and AC 


AB * m&M 

COS (p 


(4.12a) 


or 


r = AC =s Bcos (y-4>) 
o 2cos<j) cos ip 


Also from the same figure: 


D » AE*sin3 


and from the property of a logarithmic spiral: 


0 tan<J> 
AE = AC*e o 


(4.12b) 


(4.13a) 


(4.13b) 


By substituting Equations 4.12b and 4.13b into Equation 4.13a: 


_ Bcos(y-4>) sing 0 tan<f> 
p 2cos<j) cos\p e ° 


(4.14a) 


°3 , cos(y-<j>) sxn3 0 tan$ /, , 

= m > = ; — e o (4.14b) 

B 2cos(f> cosip 

The value of 3 is known if m' _< m (m=D/B) and is equal to the angle of 
internal friction. For m' > m, 3 will be smaller than <j). The following 
method is used to determine whether or not m' is smaller than m. 

(1) 3 » <j) is assumed. 

(2) This value is substituted into Equation 4.14b and the correspond- 
ing m' value is calculated. 

(3) If m' < m (m=D/B) , the correct value of 3 is <f>. If m 1 > m, the 
correct value of 3 is smaller than <j> (see Fig. 4.3) and can be 
found by iterative procedures. 

For 3=3 (3 =the correct value of 3 for the given variables), m' ■> m, 

c c 

or from Equation 4.14b: 


e tan+ 

^ j e o 

2cos(p*cos$ 


(4.15) 
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For given values of <j>, 6/<j), a, and D/B, the angle 3 is the only unknown in 
Equation 4.15. Thus, the following iterative procedure can be used: 

(1) From Fig. 4.3, rv= 180°-(y+if0 and AD=AC*e r|tan ^ and by substituting 
Equation 4.12b for the value of AC: 


Bcos(Y-j>) ritan<|> 
2cos(f> cosiJj 


(4.16) 


and as a first guess: 


or 


tanB = DD'/AD 
o 

3 = tan -1 

o 


2mcos4)*cosl|j 
cos (y- 4>) *e^ tan< ^ 


(4.17) 


Thus, 3 can be calculated explicitly from Equation 4.17 because n is known, 
o 

(2) This value of 3 = 3 is substituted into: 

o 


sinB 

n 


2mcos(l>*cosi|J 
co s ( y-<}> ) e^o tan< ^ 


where 6 = 180° - (y+^)+3 and the new value of 3, (3=3 ) is 

o o n 

calculated from: 


3 


n 


sin -1 


2mcos(j)*cos^ 
cos (y-<t>)e®o tan< ^ 


(4.18) 


(3) Iteration is stopped after the first step if: 

| 3 - 3 | < 0 . 1 ° 

.. 'no 1 — 

and the value of 3 is taken equal to 3 • If not, for the second 
n c 

iteration; 


3 


o 


B o + 3 
2 


n 
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arithmetic average of first guess and first iteration values is 
used, and the same procedure is repeated until the error in 3 q 
is equal or smaller than 0.1°, 

In most cases only few iterations were necessary. Fig. 4.4 is an example of 
the convergence of this iterative procedure. 


Procedure 

Once 3 is known, the factors N and N can be calculated from Equations 

c Yq 

4.6 and 4.9, respectively, for given values of <f), a, 6/<j>, and D/B. A 
computer program was written to calculate and N separately for different 
values of these four parameters. The complete program consist of Program 
NC Program NGQ and one subroutine, ANG which calculates angle y (the topmost 
angle of the plane shear zone) . The angle 3 Q is calculated in the main 
program. A flow diagram for the program is given in Fig. 4.5, and a listing 
is given in Appendix A. 

Example curves showing and N versus 4> for a=15° to 90°, 6/(j>-0.0, 

0.5, and 1.0, K=K o =l-sin<j), and for D/B=10.0 are given in Figs. 4.§ and 4.7. It 
should be noted that the break points in Fig. 4.6 for are a function of 

relative depth. The dashed lines indicate the versus (J) relationship for 

higher D/B values than indicated on the figure. For the relative depths 
indicated, solid lines should be used. A complete set of curves for several 
values of a, 6/(j), and D/B for general shear failure is presented in Appendix 
A. 


Effect of Base Apex Angle 

Theoretical and experimental (see Chapter Six) results both show that 
the ultimate penetration resistance of soil to wedge (or cone) shaped penetro- 
meters depends strongly on the base apex angle (2a). Fig. 4.8 shows the 
effect of the base semiapex angle on the bearing capacity factors and 
N for various base roughnesses, for given values of soil friction angle 
(cj)=30 o ) and relative depth (D/B=10.0). 

For rough bases, Fig. 4.8 indicates that the bearing capacity factors 
do not change for semiapex angles larger than approximately 15°. This is due 
to the fact that a soil wedge (or cone) develops in front of blunt rough 
bases during penetration. However, the bearing capacity factors for rough 
bases do increase with decreasing values of a below 15°. 

For smooth bases, the bearing capacity factors are strongly dependent 
on the value of the base semiapex angle. Fig. 4.8 indicates that the values 
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Angle, jS -degrees 



FIG. 4.4 AN EXAMPLE OF THE CONVERGENCE 
OF ITERATION PROCEDURE FOR 
THE CALCULATION OF ANGLE £ c 
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FIG. 4.6 Nc FACTORS FOR D/B = 10.0 
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FIG. 4.8 VARIATION OF BEARING CAPACITY FACTORS WITH BASE SEMIAPEX ANGLE 




of the factors increase sharply with increasing apex angles*. 

Effect of Base Roughness ■ . 

The bearing capacity factors N £ and H are dependent not only on a, 
but also on the base roughness (6/<j>). Therefore, knowledge of a alone is 
insufficient information from which to draw conclusions regarding penetra- 
tion resistance. 

As presented in Chapter Two, the base roughness has little influence 
on the penetration resistance encountered by a plane and horizontal contact 
area in cohesive (<f>=0) soils. However, at the surface of a cohesionless 
soil, a smooth flat-ended penetrometer will encounter less penetration 
resistance than one having a rough base. The effect of roughness on the 
penetration resistance of flat-ended fjenetrometers in cohesionless soils 
diminishes for relative depths greater than 4.0 as can be seen from Figs. 

A. 7 through A. 10. 

It can be seen from Fig. 4.8 that the effect of base roughness is 
specially important for sharp (small a) wedge or cone shaped penetrometers. 
However, this effect vanishes for larger values of base semiapex angle. 

Fig. 4.9 shows the variation of the bearing capacity factors N c and 
N with base roughness for a base semiapex angle of 15° . It can be seen 

i HI 

that the bearing capacity factors for a given intermediate roughness should 

not be estimated by linear interpolation between perfectly smooth and 

perfectly rough values. However, linear interpolation between 6/<j>= 0.0 

(perfectly smooth) and 6/<p= 0.5 (semi-rough) and between 6/<p= 0.5 and 

6/<j>= 1.0 (perfectly rough) yields a very good approximation to the correct 

values of N and N 

c yq 

Effect of Initial Stresses 

It can be seen from Equation 4.9 that the bearing capacity factor N 
is dependent upon the value of lateral earth pressure coefficient (K) . The 
value of K is determined by the penetrometer shape (relative sizes of 
penetrometer base and shaft), soil density, and initial stresses. An increase 
in the initial lateral stresses means an increase in the value of K, It has 
been shown that, for example, vibration of cohesionless soils will cause an 
increase in the value of K (D 'Appolonia, et al. , 1969). Therefore, the 
effect of initial lateral stresses on penetration resistance can implicitly 
be studied by assigning different values for K. As an example. Fig. 4.10 
shows how the value of K increases with increasing values of K. 
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Roughness of penetrometer base, Roughness of penetrometer base, 

FIG. 4.9 VARIATION OF BEARING CAPACITY FACTORS WITH BASE ROUGHNESS 
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FIG. 4.10 EFFECT OF INITIAL STRESSES ON PENETRATION RESISTANCE 








Comparison with Other Available Solutions 

As discussed in Chapter Two, Brinch Hansen (1961) proposed that bearing 
capacity factors corresponding to D/B=0 be modified by depth factors in 
order to calculate bearing capacities corresponding to larger relative depths. 
Based on that study and investigations by Meyerhof (1951, 1955), Skempton (1951) 
and others, Brinch Hansen proposed the following generalized, semi-empirical 
depth factors: 


d 

c 


1.0 + 


0.35 


h 


0.6 


(l+7tan 4 <f>) 


d 

q 


d 


c 


d -1.0 

c 


N 

q 


d = 1.0 
Y 


(4.19a) 


(4.19b) 


(4.19c) 


It is obvious that d c and d^ approach limiting values with increasing 
relative depth for a given value of <f>. Also, for <j)>25 0 , Brinch Hansen (1961) 
noted that the second term of Equation 4.19b approaches zero. Thus, for 

<j»25°: 


d = d (4.20) 

9 e 

For high values of relative depth, the factor d £ approaches an ultimate value 
given by: 


(d c^ult. = 1,0 + 0.6/(l+7tan 4 <|>) 


(4.21) 


Ultimate, depth factors piay also be calculated theoretically from 
Meyerhof’s (1961a) values of N c corresponding to shallow and deep founda- 
tions by taking the following ratio: 


N ("deep" foundations) 

\ “ £ , 

'c'ult. N c ("shallow" foundations) 


(4.22) 
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Further, the theory proposed in the current investigation may also be used 
to formulate ultimate depth factors as follows: 

N (for D/B>(D/B) ) 

^c^ult. = N c (for D/B = 0) (4.23) 

It should be noted that in Equation 4.23 the value of remains constant 

for relative depths equal to or greater than the critical relative depth. Cal- 

} 

culated values of the ultimate depth factors for the cohesion term obtained 
from Equations 4.21, 4.22, and 4.23 as a function of <j> are compared in 
Fig. 4.11. It can be seen from Fig. 4.11 that there is a very close agree- 
ment between the values calculated from the theory developed in this 
investigation and the values proposed empirically by Brinch Hansen (1961) . 

The influence of base roughness and base apex angle on the ultimate 
depth factors were also investigated. It can be seen from Fig. 4.12 that 
the effect of roughness is relatively insignificant. 

BEARING CAPACITY OF CIRCULAR PENETROMETERS 

The calculated bearing capacities are applicable to strip foundations. 

In practice, however, most penetrometers are circular in shape. Further, as shown 
in Chapter Three, the deformations around a circular penetrometer are 
different than those observed for strip foundations. The failed soil zone 
was restricted to a smaller volume and no distinct failure surfaces were 
observed in the case of the circular penetrometers. 

There have been some attempts (Larkin, 1968, Nowatzki, 1971, and 
Nowatzki and Karafiath, 1972) to treat this three-dimensional problem the- 
oretically by assuming that the Haar and von Karman hypothesis is valid; 
i.e., by assuming that the cirqumf erential stress is equal to the minor 
principal stress. However, according to Hansen and Christensen (1969), 
theoretical calculations of the bearing capacities of circular footings or 
penetrometers are open to doubt because the Haar and von Karman hypothesis 
is incompatible with the proportionality between the plastic strain and 
deviator stress tensors. Also, the assumed slip surface; i.e. , convergence 
back upon the penetrometer shaft (Nowatzki, 1971) with a discontinity at 
base level (Nowatzki and Karafiath, 1972) is incompatible with that observed 
in model tests. 
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FIG. 4.11 COMPARISON OF ULTIMATE DEPTH 
FACTORS FOR COHESION TERM 


88 




FIG. 4.12 INFLUENCE OF ROUGHNESS OF BASE ON 
ULTIMATE VALUES OF DEPTH FACTORS 
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For these reasons it was concluded that theoretical calculations of the 
bearing capacities of circular foundations are still so uncertain that best 
estimates are made by the application of shape factors, found empirically, 
to the bearing capacity factors for strip foundations. 

Many empirical shape factors have been proposed. Those recommended by 
Brinch Hansen (1961) were adopted in this study, because they agreed very 
closely with the experimentally determined values (see Chapter Six) , The 
following equations, determined semi-empirically, were presented by Brinch 
Hansen (1961) for the individual shape factors: 


E c - 1.0 + (0.2 + tan 6 <(>) (B/L) 

(4.24a) 

5 c " 1 


E - E - c 
’q S N 

(4.24b) 

q 


Ey = 1.0 - | (0.2 + tan 6 <}>) (B/L) 

(4.24c) 


Brinch Hansen (1966) later suggested that Ey be given by: 

E ~ 1.0 - 0.4 (B/L) (4.25) 

Brinch Hansen also pointed out that for <j) > 25° : 


E - E - 1.0 + [ (0.2 + tan 6 <j>) (B/L) ] 
q c 


(4.26) 


It may be seen from Equation 4.26 that the shape factor is inde- 
pendent of relative depth. However, because the contributions of friction 
and surcharge were considered together in formulating the bearing capacity 
factor N , it may be shown that the corresponding shape factor § is 
dependent on relative depth. 

The procedure discussed below is used for the calculation of the shape 


factor E. 


yq 


As previously discussed, the unit resistance of strip and 


circular (or other than strip in shape) foundations may be expressed as: 


(q f ) strip 


1 

2 


y BN 
's y 


+ Y DN 
s q 


(4.27a) 
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(4.27b) 


or 


(q J , - 4 Y BN £ + Y DN £ 

VH f' general 2 ' s y Y s q 9 


strip -f\ + t H q 


(4*28a) 


qf / general 


1 N £ + ^ N £ 

2 Y^y B q s q 


(4.28b) 


The ratio of q f / (Y 8 B) general to q f / (Y g B) strip is the shape factor g yq and 
is given by: 



4 N £ + mN £ 

2 h q q 

N + mN 

2 Y q 


(4.29) 


where m = D/B. 

At the surface (D/B=0) and the contribution of surcharge due to over- 
burden will be zero, consequently; 

£ - £ - 1.0 - 0.4 (B/L) (4.30) 

yq Y 


or 



= 0.6 for circular penetrometers. 


For large values of relative depth, however, the contribution of the 
friction term can be neglected compared to the contribution of the surcharge 
(overburden) term, yielding the following expression from Equations 4.26 and 
4.29: 


£ = £ - 1.0 + (0.2 + tan 6 !})) (B/L) (4.31) 

yq q 

Values of the shape factor .£ may be calculated from Equation 4.29 for 
different values of <j> and relative depth (m=D/B) . The factors £^ and £ q are 
determined from Equations 4.25 and 4.26, and the bearing capacity factors N^ 
and N are taken from Fig. 2.4. The calculated values of £ are presented 

q yq 

in Fig. 4.13. 

However, in order to apply the procedure suggested in Chapter Seven 
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Relative depth, % 



FIG. 4.13 VARIATION OF SHAPE FACTOR £yq WITH DEPTH 
FOR GIVEN VALUES OF ANGLE OF INTERNAL 
FRICTION 


92 



’yq 


for the determination of in-situ strength properties from the results of 
cone penetration tests, it is necessary to formulate the variation in 
analytically. The following empirical relationship has been obtained which 
agrees closely with the values calculated theoretically from Equation 4.29 
for the range of interest: 



= 1.0 - 0.4(B/L) + 



1.5 

1.5 


(0.6+tan 6 <J>) (B/L) 


(4.32) 


where 


£ yq - shape factor for the friction-surcharge term 
B = width of penetrometer base 

D = depth of penetrometer base 

L = length of penetrometer base ( -B for circular penetrometers) 
<f> = angle of internal friction of soil 

At the surface, (D/B=0) , Equation 4.32 becomes: 


Sy *■ ? Y - 1.0 - 0.4 (B/L) (4.33) 

Further, it can be seen that, for large values of relative depth (D/B), 

Equation 4.32 becomes: 

S yq = £ = 1.0 + (0.2 + tan 6 <f>) (B/L) (4.34) 

Therefore, it can be seen that Equation 4.32 satisfies known relationships 
(Equations 4.25 and 4.26) at the limits (D/B = 0 and *>) . 

Table 4.1 presents a comparison (over a wide range of values of <j> 
and D/B) of values of 5 yq by Equation 4.32 and by Equation 4.29. It can be 
seen that the agreement is quite good, and thus it may be concluded that 
Equation 4.32 may be satisfactorily utilized in the analysis of cone penetra- 
tion data. 


SUMMARY 

Knowledge of the failure mechanism associated with static penetration 
(see Chapter Three) has permitted development of a new relationships for 
ultimate resistance to penetration of wedges which account explicitly for base apex 
angle, base roughness, soil friction angle, and relative depth. The derivation of 
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Note: Upper values are by Equation 4.29, lower values are by Equation 4.32. 


















the bearing capacity factors and N in terms of these parameters and 
curves relating these factors to the angle of internal friction of soil for 
several values of a, S/(j>, and D/B for strip foundations were presented. 

The values of ultimate depth factors calculated from these relationships 
agree closely with the empirical values given by Brinch Hansen (1961) . The 


shape factor 
of E and £ 

Y ■ q 



was calculated in terms of relative depth using the values 


proposed by Brinch Hansen (1961) . Based on the results of 


these calculations, an empirical relationship was proposed which relates 


§ to the angle of internal friction of soil and relative depth. 


95 



CHAPTER FIVE 


EXPERIMENTAL INVESTIGATION - SOIL PROPERTIES 


INTRODUCTION 

In order to evaluate the proposed theory experimentally, accurate 
values of the angle of internal friction of the soil and the penetro- 
meter to soil friction must be determined. Comprehensive series of drained 
triaxial and direct; shear tests were performed to measure the strength 
properties of the test soil, and a series of direct shear tests was done 
in order to determine the penetrometer to soil friction. The results of 
these tests enabled the formulation of: 

(1) The soil friction angle (<j>) in terms of the pre-shear void ratio 
and mean normal stress (or confining pressure) , 

(2) The penetrometer to soil friction angle (6) in terms of the pre- 
shear void ratio and normal stress for various penetrometer 
materials and, 

(3) The variation of roughness (<$/<))) with pre-shear void ratio. 

These results were used to predict the penetration resistance of the 

soil using the proposed theory, thus enabling a comparison between predicted 
and measured values (see Chapter Six) . 

SOIL CLASSIFICATION 

A fine, clean, air-dried sand (Monterey Sand No. 0) was chosen for 
this study. Classification data for this sand are summarized below: 


Mean diameter 0.36 mm 

Coefficient of uniformity 1.45 

Specific gravity of grains 2.648 g/cm 3 

Maximum void ratio 0.825 

Minimum void ratio 0.558 


The gradation curve of Monterey Sand No. 0 is given in Fig. 5.1. Microscopic 
observation of the sand showed that the particles were sub-rounded and that 
quartz and feldspar were the predominant minerals with a trace of mica. 

STRENGTH PROPERTIES FROM TRIAXIAL COMPRESSION TESTS 
Test Procedure 

Specimens were prepared in a triaxial cell to the desired uniform 
initial densities. The desired chamber pressure (pressures used were 
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FIG. 5.1 GRADATION CURVE OF MONTEREY SAND No. 



29.4, 58.8, 98, and 294 kN/m 2 ) was applied by air pressure on the chamber 
water and the specimen was tested 30 minutes later in a Wykeham-Farrance 
strain-controlled testing machine at a deformation rate of 0.006 inches 
per minute (0.152 mm/min) . During the test, measurements were taken of 
the axial load with a load cell, axial deformation with a dial gage and 
change in volume of chamber water with a volume change device of the type 
developed by Chan and Duncan (1967). 

Test Results 

A total of 16 drained triaxial shear tests was conducted. The results 
of these tests are plotted in Figs. 5.2 and 5.3 which show the values of 
principal stress ratio (a^/a^) versus axial strain for the conditions 
investigated. The peak values of soil friction angles were calculated 
from the results of these tests using the equation: 

<t> = sin" 1 {((a 1 /a 3 ) f - U/[(a 1 /a 3 ) f + 1]} (5.1) 

The triaxial shear tests were carried to axial strain values of 10 
to 15 percent. In the determination of residual triaxial friction angles, 
principal stress ratios corresponding to axial strains of from 8 to 12 
percent were used. The change in principal stress ratio beyond these 
axial strain values was negligible for the range of confining pressures 
and soil densities used. 

Peak and ultimate triaxial friction angles are presented in Table 
5.1 and are plotted against pre-shear void ratio (e £ ) for the various 
confining pressures in Fig. 5.4. Values of the peak friction angles deter- 
mined from plane strain tests by Lade (1972) for the same sand are also 
plotted on Fig. 5.4 for comparison. 

It can be seen that the difference between peak and ultimate values 
of friction angles obtained from triaxial shear tests is not great, with 
almost no difference for loose soil conditions. 

Peak principal stress ratios were found to decrease with increasing 
confining pressures (or mean normal stress) for a given pre-shear void 
ratio. As a result, the failure envelopes shown in Fig. 5.5 are slightly 
curved, indicating that an increase in mean normal stress causes a reduc- 
tion in soil friction angle values. 
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Axial strain, € Q -(%) 

FIG. 5*2 STRESS -STRAIN RELATIONS! 
DETERMINED FROM TR1AXI; 
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Axial strain, € 0 -(%) Axial strain, € a -(7< 

FIG. 5.3 STRESS-STRAIN RELATIONSHIPS OF AIR - DRY MONTEREY SAND No. 
DETERMINED FROM TRIAXIAL SHEAR TESTS 





Angle of internal friction, <£- degrees 
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0.50 0.55 * 0.60 0.65 0.70 0.75 0.80 0.85 0.90 

Void ratio, e 

FIG. 5.4 SUMMARY OF FRICTION ANGLES FOR MONTEREY 
SAND No. 0 
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Table 5.1 SUMMARY OF TRIAXIAL PEAK AND ULTIMATE FRICTION 
ANGLES FOR MONTEREY SAND NO. 0 


Confining 

Pre-shear 

Peak Friction 

Ultimate 

Pressure 

Void Ratio 

Angle 

Friction 

C^-kN/m 2 

e 

<f> -deg, 

Angle. 

o 

c 

P 

4> r -deg. 


0.802 

36.4 

35.9 


0.760 

39.0 

— 

29.4 

0.692 

41.2 

39.8 


0.585 

46.5 

44.0 


0.802 

36.0 

35.7 


0.736 

39.0 

— 

58.8 

0.692 

40.9 

- 


0.585 

46.5 

44.5 


0.797 

35.1 

35.0 


0.735 

38.0 

_ 

98.0 

0.669 

41.2 

- 


0.585 

45.2 

42.3 


0.783 

33.6 

33.5 


0.731 

35.8 

— 

294.0 

0.681 

37.3 

36.5 


0.575 

42.1 

40.6 
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FIG. 5.5 FAILURE ENVELOPES FOR MONTEREY SAND No. 0 




As suggested by Banks (1968) , the triaxial <f> values were plotted 
against the logarithm of confining pressure as shown in Fig. 5.6. It may 
be seen that the relationships can be represented by straight lines for 
the range of confining pressures used. This variation of <j) with may be 
expressed in the following form: 

<K e c » a 3 ) “ -m(e c ) log (a 3 /p a ) (5.2) 

in which ^ (e^) is the angle of internal friction at a given pre-shear void 
ratio and at a confining pressure of p (atmospheric pressure) expressed 
in the same units as cr_, and m(e ) is the reduction in <J> per log cycle 
change in confining pressure. For the range of confining pressures 
(29.4 to 294 kN/m 2 ) used, m values are between 3.5 and 6.0. The triaxial 
friction angle of Monterey Sand No. 0 can be determined from the Figs. 5.4 
and 5.6 for any value of confining pressure (or normal stress on the 
failure surface) and pre-shear void ratio. 

STRENGTH PROPERTIES FROM DIRECT SHEAR TESTS 

Penetrometer to soil friction angles (<5) were determined from the 
results of direct shear tests. A series of direct shear tests was 
also done on the sand in order to obtain values of the soil friction 
angle for determination of 6/<J>, 

The results of soil to soil direct shear tests are shown in Fig. 

5.7 and indicate the variation of shear stress with relative shear 
displacements for normal stresses of 168, 353, and 535 kN/m 2 . The 
corresponding direct shear friction angles were calcualted and are pre- 
sented in Table 5,2. These values can be compared with the values obtained 
from triaxial tests, as shown in Fig. 5.8. Because the values are stress 
dependent, comparisons should be made at same normal stress values. The 
normal stress on the failure plane in a triaxial test is calculated from 
the following relationship: 

O mf/ 0 3f = 1 + sin ^ (5.3) 

Fig. 5.8 shows the average <j> versus void ratio relationship which will be 
used for the calculation of 6/<j> in the following section. 
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FIG. 5.6 VARIATION OF TRIAXIAL FRICTION ANGLE OF MONTEREY SAND No. 
WITH CONFINING PRESSURE 
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FIG. 5.7 SHEAR STRESS VERSUS SHEAR DISPLACEMENT 
RELATIONSHIPS OF MONTEREY SAND No. 0 
DETERMINED FROM DIRECT SHEAR TESTS 





FIG. 5.8a COMPARISON OF SOIL TO SOIL TRIAXIAL AND 
DIRECT SHEAR FRICTION ANGLES 



FIG. 5.8b SOIL TO SOIL FRICTION ANGLES OF MONTEREY 
SAND No. 0 FROM DIRECT SHEAR TESTS 
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Table 5.2 SUMMARY OF SOIL TO SOIL DIRECT 
SHEAR TESTS 


Normal Stress 
on failure plane 
a mf -kN/ m 2 

Initial Void 
Ratio , e 

Angle of 
Internal 
Friction 
<|)-deg. 


0.748 

35.0 

168 

0.632 

42.6 


0.550 

50.1 


0.759 

34.2 

353 

0.608 

43.2 


0.555 

47.8 


0.761 

34.8 

535 

0.619 

41.4 


0.560 

46.0 

— — 


PENETROMETER TO SOIL FRICTION 

To determine the friction angle between the penetrometer material and 
soil, and establish a basis for the variation of 6/<j> with void ratio, tests 
were conducted in a similar manner to the standard direct shear test. The 
upper half of the shear box was filled with Monterey Sand No. 0 at a desired 
initial density, and the lower half of the shear box was replaced by a solid 
sample of penetrometer material. A section through the shear box with the 
samples in place is shown in Fig. 5.9. The three differenct penetrometer 
materials tested are listed below: 

(1) Polished hard-anodized aluminum 

(2) Hard- anodized aluminum 

(3) Sanded aluminum (sand particles glued to aluminum) . 

The tests were conducted in a Karol-Warner direct shear machine connected 
to a variable speed motor in series with a gear reduction box so that the rate 
of shearing could be maintained at 0.0015 in/min (0.048 mm/min). A total of 
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FIG. 5.9 DIRECT SHEAR TEST SETUP 



15 tests at various void ratios was conducted using normal pressures of 168, 

353, and 535 kN/m 2 , and the results are summarized in Table 5.3 and presented 
graphically in Fig. 5.10. The following conclusions may be drawn based on 
the test results: 

(1) A wide range of penetrometer to soil friction angles may be obtained 
depending on the penetrometer material. 

(2) Penetrometer to soil friction angles are strongly affected by the 
initial void ratio of the soil, as are also the spil to soil friction 
angles . 

(3) Penetrometer to soil friction angles decrease slightly with increas- 
ing values of normal stress on the failure surface as ih the case of 
soil to soil friction tests. However, this effect can be neglected 
by employing an average value of (6) for the range of normal stress 
values . 

The friction angle 6 between various penetrometer materials and the soil, 
as well as the angle of internal friction (<|)) of the soil can be determined 
from Figs. 5.8 and 5.10 for given values of soil void ratio. Values of 6 / <J> 
were calculated and are presented in Table 5.4 and Fig. 5.11. It can be seen 
that an essentially unique value of 6/<p exists regardless of the soil density 
or soil friction angle (<(>). Values of <5/<j> were found to be approximately 
equal to 0.3, 0.5, and 0.9, respectively, for polished hard-anodized aluminum, 
hard-anodized aluminum, and sanded aluminum. 

SUMMARY 

The internal friction angle of Monterey Sand No. 0 was determined from 
the results of triaxial and direct shear tests. In addition, a series of direct 
shear interface tests has been performed to determine the penetrometer to soil 
friction for three different penetrometer materials. Using these results, 
both 4> and 6 have been formulated in terms of soil void ratio over a range of 
normal stresses. These results are used in the next chapter to predict the 
penetration resistance of Monterey Sand No. 0 in accordance with the relation- 
ships developed in Chapter Four. 
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Table 5.3 SUMMARY OF PENETROMETER TO SOIL FRICTION 
TESTS IN DIRECT SHEAR 


Nature of 
Friction 

Normal Stress 
on failure plane 
CT^f-kN/m 2 

Initial 
Void Ratio 
e 

Friction 

Angle 

6-deg. 



0.776 

10.1 


168 

0.680 

11.9 



0.510 

18.4 



0.722 

9.7 

Soil to 




polished 

353 

0.665 

10.8 

hard- anodi z ed 




aluminum 


0.551 

15.4 



0.792 

9.1 


535 

0.672 

10.5 



0.545 

13.2 



0.760 

16.1 

Soil to 




hard-- anodized 

353 

0.620 

20.9 

aluminum 






0.555 

23.8 



0.743 

31.7 

Soil to 




sanded 

353 

0.587 

40.2 

aluminum 






0.535 

44.2 
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FIG. 5.10 PENETROMETER TO SOIL FRICTION ANGLES 
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Table 5.4 SUMMARY OF ROUGHNESS (S/<|>) VALUES FOR DIFFERENT PENETROMETER MATERIALS 
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Penetrometer to soil friction 




CHAPTER SIX 


EXPERIMENTAL INVESTIGATION - 
LABORATORY PENETRATION TESTS 


COHESIONLESS SOIL 
Introduction 

In order to check the accuracy of the theoretical predictions, carefully 
controlled laboratory tests were conducted. This section presents the results 
of controlled penetration tests using Monterey Sand No. 0 and comparisons of 
these results with predicted values. 

Equipment, Facilities, and Penetrometers 

Penetration tests were conducted in a test box 3.5 feet (-1.05 m) wide, 

7.0 feet (~ 2 . 1 m) long, and 4.0 feet (-1.2 m) deep. The wooden test box is made 
up of two 2-foot high removable sections. The loading equipment included a 
hydraulically activated loading apparatus consisting of a 4.0 inch (102 mm) 
diameter double acting cylinder and a flow control system. The load was trans- 
mitted to the penetrometer by a shaft rigidly connected to the piston. The 
system was designed for pressures up to 150 psi (-1035 kN/m 2 ), and has a 
maximum loading capacity of 1800 pounds (-7650 N) . The available piston stroke 
was 20.0 inches (0.508 m) . (For deep penetration tests, a different cylinder 
with a stroke of 5 feet (-1.50 m) was used). The hydraulic system allowed 
application of the load either incrementally or continuously at any desired 
loading rate or at a constant penetration rate. 

The penetration test equipment also included a supporting frame consisting 
of two aluminum channels spanning across the test box. The heavy supporting 
frame acted as a reaction for the loading piston. The supporting beam for 
the cylinder is connected by bolts to two channeled aluminum columns, which 
are individually supported on wheels and are therefore easily movable. The 
elevation of the supporting beam is easily adjusted to allow testing of any 
depth of soil deposit in the test box. 

The two basic types of penetrometer tips used in the laboratory penetra- 
tion tests were rectangular (strip) and circular in section. The rectangular 
tips had a length-to-width ratio of approximately 6 to 1. The tips were 
detachable from the shaft and had various semiapex angles and roughness. A 
complete listing of the penetrometer tips used in the tests is given in Table 6.1. 
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Table 6.1 SUMMARY OF PENETROMETER TIPS USED IN PENETRATION 
TESTS OF MONTEREY SAND NO. 0 


Section 

Shape 

Base semi-apex 
Angle, a 

Tip 

Material 

Roughness 

(6/<J>) 


15 

Hard-anodized 

aluminum 

0.5 


Sanded 

aluminum 

0.9 

Rectangular* 

30 

Hard-anodized 

aluminum 

0.5 


Sanded 

aluminum 

0.9 


90 

Hard-anodized 

aluminum 

0.5 



Polished hard 
anodized 
aluminum 

0.3 


15 

Hard-anodized 

aluminum 

0.5 



Sanded 

aluminum 

0.9 

Circular** 

30 

Polished hard 
anodized 
aluminum 

0.3 


Sanded 

aluminum 

0.9 


90 

Hard-anodized 

aluminum 

0.5 


Sanded 

aluminum 

0.9 


*Tips have dimensions of 0.8" x 5.0" (20.3 mm x 127 mm), and penetrometer 
shaft dimensions are 0.625" x 5.0" (15.9 mm x 127 mm). 

**Tips are 0.8" (20.3 mm) in base diameter, and the shaft is 0.625" (15.9 mm) 
in diameter. 
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Samples of different but uniform densities which were reproducible from 
test to test were required. While many sample preparation techniques are 
available (screen technique for loose, vibration technique for dense samples), 
pluvial compaction offers one of the most convenient means of preparing large 
uniform sand deposits at any required relative density. In pluvial compaction, 
air-dry sand is allowed to fall through the air in order to build up the 
required layer. Experiments by Kolbuszewski (1948) have shown that the factors 
controlling the final density are height of fall and intensity of sand rain. 

The method has some limitations , because it tends to produce layering in well- 
graded materials and it induces some cross anisotropy. The problem of layering 
can be eliminated by using uniformly graded materials, (The uniformity co- 
efficient of the sand tested is around 1.5), 

In addition to the ability of pluvial compaction to deposit uniform sand 
layers, the method also allows the deposition of large quantities of sand 
quickly and easily. In fact, the entire operation can be automated with a sand 
spreader box such as that shown in Fig. 6.1. The spreader box is similar to 
the device described by Walker and Whitaker (1967) and used by Silver (1970), 
and consists of a 1/3 cubic yard (0.26 m 3 ) wooden hopper sloping down to an 
opening at the bottom which is closed off by an aluminum roller. A rubber 
seal prevents sand from falling from the front of the box while the roller is 
stationary. The roller is driven through a chain drive by a variable speed 
motor. Limit switches make it possible to start and stop the roller at pre- 
determined locations. 

The spreader box is mounted on a frame that spans over the test box. The 
frame is chain-driven by an electric motor. 

During operation, the roller rotates and sand is ejected from the box in 
a uniform sheet as shown in Fig. 6.1. The rate of deposition can be varied 
by changing the roller speed. The uniform rain of sand, the constant forward 
speed of frame, and the constant height drop cause uniform layers of sand to 
be deposited. As the soil deposit builds up, the spreader box is elevated to 
maintain a constant height drop. 

Measurements of sand layer density showed that relative density did not 
vary more than 5 percent from the desired value throughout the layer (see 
Appendix B) . For the sand used in this investigation, it was found that rela- 
tive densities of from 35 to 100 percent could easily be obtained using the 
spreader box. 
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FIG. 6.1 SPREADER BOX AND 
SAND DEPOSITION 
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Test Procedure 


Soil densities ranged between 35 to 100 percent relative density. The 
test box was large enough to allow 4 to 5 separate penetration tests on one 
sample. Penetrations were spaced so as to minimize test box side effects and 
the influence of other penetration tests. The load transfered to the penetro- 
meter was measured with a load cell and recorded with an electronic recorder- 
print out system. The skin friction resistance of soil along the penetrometer 
shaft was insignificant due to smaller shaft diameter than the base diameter. 
The vertical penetration was measured by observing the relative movement of a 
fixed marker on a ruler. 


Prediction Method 
General 

As described in Chapter Four, the unit penetration resistance can be 
determined for cohesionless soils by: 


q £ = Y BN £ 

Vs Tq Yq 

in which q^ = unit penetration resistance 
y = unit soil weight 


N = bearing capacity friction- surcharge factor 

£ = shape factor 

yq F 

B = width of penetrometer base 


( 6 . 1 ) 


In dimensionless form, Equation 6.1 may be written; 


q f /(Y s B) = 


N 5 

yq yq 


( 6 . 2 ) 


Equation 6.2 indicates that, because both and are functions of relative 

depth, the penetration data should be presented in the form of q f /(Y s ) versus 

x s 

D/B. 


Estimation of mean normal stress 

As discussed in Chapter Two, the approximate average value of mean normal 
stress along a failure surface can be estimated from: 


°g,M ' V 10 


(6.3) 


in which a 

g,M 

resistance. 


is the average mean normal stress. 


and is the unit penetration 
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Determination of proper soil friction angle 

The penetration problem in cohesionless soils is complicated by the 
fact that during the gradual load increase on the soil the shear strength is 
not immediately mobilized at all points of the slip surface, but at first only 
at the points where the shearing stresses are largest, with gradual extension 
to other points. This gradual progression causes modifications of the soil 
properties along the slip surface, especially in dense soils. In highly 
stressed zones, the soil begins to dilate and a decrease of density and thus 
a decrease in shear strength takes place. Therefore, when the state of 
rupture along a slip surface is reached, the shear strength corresponding to 
the original density is not available along the whole surface. 

Further, the stress-strain characteristics of the soil must be considered 
in addition to the strength properties. Typical stress-strain relationships 
for plane strain and triaxial tests and dense and loose sands are shown in 
Fig. 6.2. It can be seen that if deformations at some points along the failure 
surface before rupture are sufficiently large to surpass the peak points of 
the stress-strain curves, the shear strength will drop to ultimate (residual) 
values. Consideration of stress-strain characteristics is most important in 
the case of a dense sand under plane strain conditions, because the stress- 
strain curve has a very distinct peak as shown in Fig. 6.2. Consequently, 
the use of friction angles corresponding to plane strain peak values should 
cause overestimation of the penetration resistance of cohesionless soils. 

Because the ultimate values of shear stresses control the resistance at 
large deformations, friction angles corresponding to ultimate values of shear 
stresses (ultimate friction angle) should be used in the analysis of the 
continuous penetration problem. 

It is common practice in shear testing to obtain only peak values of the 
friction angle and to terminate the test before accurate ultimate values are 
obtained. However, as shown in Table 5.1, the triaxial peak and ultimate, values 
for loose sands are the same and, for dense sand, the ultimate values are only 
slightly less than peak values. Therefore, as shown in a subsequent section, 
the use of peak triaxial friction angles for dense sands should cause only a 
small overestimation of the unit penetration resistance. 

Procedure 

The penetration resistance was predicted for different values of relative 
depth using the following procedure. For a given test, the following information 
is known: 
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Principal stress ratio 




FIG, 6.2 TYPICAL STRESS -STRAIN CURVES FOR 
COHESIONLESS SOILS 
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(1) Unit weight (or void ratio) of soil deposit, 

(2) Unit penetration resistance (q^) versus penetration (D) relationship. 
The following information is needed for the calculation of penetration resistance 
using Equation 6.1: 

(1) From Equation 6.3 the average mean normal stress corresponding to a 
given value of relative depth (D/B) , 

(2) From Fig. 5.4, the angle of internal friction of the soil corresponding 
to given (determined) values of void ratio and normal stress, 

(3) From Equation 4.9, the bearing capacity (penetration resistance) 
factor N as a function of a, <p, <$/$, and D/B, 

yq 

(4) From Fig. 4.13 shape factor as a function of <f) and D/B. (For 
wedges a shape factor of 1.0 is used). 

Test Results and Predictions 

The results of penetration tests performed using wedge-tipped penetro- 
meters are shown in Figs. 6.3 through 6.7 in the form of unit penetration 
resistance versus penetration depth and also in the dimensionless form of re- 
sistance parameter q f / (y B) versus relative depth (D/B). Predicted values using 
peak values of triaxial friction angles are also shown in Figs. 6.3 through 6.7. 

The results of penetration tests performed using cone-tipped penetrometers 
are shown in Figs. 6.8 through 6.11 in the form of unit penetration resistance 
versus penetration depth and also in the form of the resistance parameter 
q^/ (Y g B) versus relative depth (D/B) . Predicted values using peak values of 
triaxial friction angles and Equation 4,32 for the shape factors are also shown 
in Figs. 6.8 through 6.11. 

It may be seen in Figs. 6.3 through 6.11 that the agreement between pre- 
dicted and measured values is quite good. This suggests that the proposed 
analytical method may be confidently used to predict the penetration resistance 
of cohesionless soils, at least to the relative depths tested. 


Determination of Shape Factors 

The bearing capacity factors dfetfaStoinecT theoretically for strip foundations 
must be modified by empirically determined shape factors when computing the unit 
penetration resistance of circular penetrometers. Empirical formulas for the 
shape factor proposed by different investigators are summarized in Table 6.2. 
It can be seen that the range of calculated values is quite large. However, it 
is possible to determine the proper shape factor for Monterey Sand No. 0 by 
comparing the resistance values for cones and wedges. The average values of tlie 
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Unit penetration resistance , q^ - kN/m 2 

FIG. 6.6 COMPARISON OF MEASURED AND PREDICTED STATIC PENETRATION CURVES FOR MONTEREY SAND No. 0 
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Table 6.2 SUMMARY OF VARIOUS EMPIRICAL FORMULAS 
FOR SHAPE FACTOR £ 

q 


Angle of 
internal 
friction 
4>-deg. 

Brinch Hansen 
(1) 

Meyerhof 

(2) 

Vesic 

(3) 

DeBeer 

(4) 

25 

1.21 

1.25 

1.67 

1.42 

30 

1.24 

1.30 

1.78 

1.50 

35 

1.32 

1.37 

1.90 

1.57 

40 

1.55 

1.46 

2.04 

1.64 

45 

2.20 

1.59 

2.20 

1.71 


(1) 

£ = 1.0 + 
q 

(0.2 + tan 6 <J>) (B/L) 

(2) 

£ - 1.0 + 
q 

0.1 tan 2 (45° + (B/L) 

(3) 

£ = 1.0 + 
q 

(0.2 + tan<j>) (B/L) 

(4) 

£ = 1.0 + 
q 

sin<j)(B/L) for D/B<1 


shape factor £ for the relative depth range (D/B=10 to 15) determined by this 
ratio procedure are shown in Fig. 6.12. The computed values of £ using 
Equation 4.32 which were derived from the equation (Equation 1 in Table 6.2) 
proposed by Brinch Hansen (1961), (see Chapter Four) are also shown in Fig. 

6.12. The close agreement between the calculated and measured values of £ 

yq 

indicates that Equation 4.32 for shape factor £ may be used for in the 
prediction of cone resistance. 

Some conclusions drawn from the test results may be summarized as follows: 

(1) The penetration resistance is sensitive to soil density. It should 
therefore be possible to use the penetration resistance of cohesion- 
less soils for estimating the in-situ soil density as well as shear 
strength properties. 

(2) Penetration resistance increases with increasing semiapex angle for 
S/(J>=0. 3. 

(3) Penetration resistance increases with increasing penetrometer roughness 

(<5/<j>). 
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Shape factor 



FIG. 6.12 COMPARISON OF MEASURED AND 
CALCULATED SHAPE FACTOR 
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(4) Penetration resistance increases with increasing relative depth (D/B), 
There was no indication of a decrease in the rate of resistance 
increase with depth even for relative depths as large as D/B=20. 

(5) Other factors being equal, the unit penetration resistance for cones 
was higher than for wedges. 

Deep Penetration Tests 

To establish the applicability of the proposed method and to investigate 
the penetration versus resistance relationships for high relative depths 
(D/B up to 40) a series of deep penetration tests was conducted. It was ob- 
served that there was no distinct change in resistance versus penetration 
relationship. The summary of the test results is given in Fig, 6.13. It can 
be seen that measured values of penetration resistance can be predicted reasonably 
well by the proposed method for relative depths as high as 40. 


COHESION-FRICTION SOIL 


Introduction 

The close agreement between predicted and measured values of the penetra- 
tion resistance of Monterey Satid No. 0 indicates the validity of Equation 4.9 
for the determination of the bearing capacity factor N for failure in general 
shear. Because lunar soil has cohesion, the validity of the proposed method 
must be further established for a cohesion-friction soil. The penetration tests 
with Lunar Soil simulant No. 2 (LSS No. 2) described in this section were 
performed in order to determine the validity of Equation 4.6 for the determin- 
ation of the value of N c » which is needed to apply the proposed method to 
such soils. At the same time the influence of soil compressibility on penetra- 
tion resistance was investigated. 


Equipment, Test Procedure, and Soil 

Static penetration tests were conducted in a test box 2.0 feet (~0.6 m) 
wide, 2.0 feet (~0.6 m) long, and 3.0 feet (-0.9 m) deep. Uniform loose and 
medium dense deposits were prepared using the constant height of drop method; 
dense deposits were prepared using a vibratory technique. The loading system 
and test procedures used were the same as previously described for cohesion- 
less soil. 

The soil used was a crushed basalt known as Lunar Soil Simulant No. 2 
having the gradational characteristics of a silty fine sand (Fig. 6.14). A 
comparison of the one-dimensional compression characteristics of this soil 
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and Monterey Sand No. 0 are presented in Fig. 6.15, which indicates that the 
simulant is much more compressible than the sand. Cohesion was obtained in 
the simulant by maintaining a small moisture content of about 2.0 percent. 

The relationship of cohesion to void ratio at various moisture contents for 
the simulant was determined by Mitchell et al. (1971) and is shown in Fig. 

6.16. 

The relationship between plane-strain ffictiop. angle of the simulant and void 
ratio is shown in Fig. 6.17 (Mitchell et al. 1971). These plane-strain values 
were converted to triaxial values according to the procedure recommended by 
Cornforth (1964). The estimated triaxial friction angle values for the simulant 
are also shown in Fig. 6.17. 

Prediction Method 

As discussed in Chapter Four, the unit penetration resistance can be ex- 
pressed for cohesion-friction soils by: 


q £ = cN £ + 

^f cc 


Y BN £ 
s Yq Yq 


(6.4) 


or, in dimensionless form: 


q f / ( Y s B ) = 


-TTT- *N £ + N £ 

Y S B c^c Yq Yq 


(6.5) 


in which q^ = unit penetration resistance 

Y = unit weight of soil 
3 

c » cohesion 
B = width of base 

N , N = bearing capacity factors f. , f«(<J>, 6/<t>, a, D/B) 
c Yq *• 

K c , 5 yq = shape factors, f 3 (<(>) , f 4 (<J> , D/B) 

The following procedure may be used to predict the ultimate penetration 
resistance of the simulant for known average values of void ratio and water 
content : 

(1) Enter Fig. 6.16 with the known values of e and w, and determine the 
cohesion value. 

(2) Enter Fig. 6.17 with the known average void ratio, and determine 
corresponding triaxial friction angle (<j>) value. 
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Compressive stress, <r - fycm 2 


FIG, 6.15 ONE- DIMENSIONAL COMPRESSION CURVES FOR LSS No. 2 
AND MONTEREY SAND No. 0 
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Cohesion , c - kN/ m z 


FIG. 6.16 RELATIONSHIP BETWEEN COHESION AND VOID RATIO 
FOR LSS No. 2 


( After Mitchell et at, 1971 ) 
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(3) 

(4) 

(5) 

( 6 ) 


(7) 


At a given relative depth (D/B) and known values of a and 6/<j) , enter 

the corresponding figures in Appendix A and determine N and N 

c Yq 

For the <j> value determined in step 2, use Equation 4,24a to 

calculate £ . 

c 

For the known values of <(> and D/B use Equation 4.32 to calcualte £ 

yq 

Substitute the values of c, N c> N * ^ c » ^yq* ^s* anc * ® into Equation 
6.4 or 6.5 to determine the ultimate penetration resistance (q^) or 
the dimensionless resistance parameter q,;/(Y B) for the particular 

r S 

D/Bo 

Repeat steps 3 through 6 for other values of D/B. 


Test Results and Predictions 

A series of penetration tests in the LSS No. 2 were conducted using a 
cone-tipped penetrometer (a=15° , B=0.8 in. , and 6/<f>=0.5). The soil was pre- 
pared over a wide range of densities in order to evaluate the applicability of 
the proposed method for both dense and loose conditions. 

Penetration test data are presented in Fig. 6.18 in the form of penetra- 
tion resistance (q^) versus relative depth (D/B) . The ratios of predicted 
to measured ultimate penetration resistance are presented in Table 6.3. It may 
be seen from Table 6.3 and Fig. 6.19 that, for dense deposits, these ratios 
are close to unity, indicating the validity of the prediction procedure and 
the values of for general shear failure Conditions. However, for low 
densities, the ratios are larger than one, indicating the significant in- 
fluence of soil compressibility on penetration resistance. This influence 
can be anticipated by noting in Fig. 6,15 that the simulant is much more 
compressible than Monterey Sand No. 0 at high void ratios. Therefore, the 
use of bearing capacity factors formulated for general shear failure condi- 
tions will cause overestimation of the penetration resistance of compressible 
deposits. As a result, of soil compressibility, the shear surface is restricted 
to a smaller zone around the penetrometer tip as shown in Fig. 6.20. Vesic 
(1963) suggested that the angle which defines the extent of the failure surface 
may be represented empirically by: 


e = r.9> 


( 6 . 6 ) 


By substituting the value of 0 determined from Equation 6.6 into Equa- 
tions 4.6 and 4.9, values of N* and N* may be obtained for compressible 
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Unit penetration resistance , - psi 



FIG. 6.18 MEASURED PENETRATION CURVES FOR LSS No.2 
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Table 6.3 RATIOS OF PREDICTED TO MEASURED PENETRATION 
RESISTANCE FOR LSS NO. 2 USING GENERAL SHEAR 
FAILURE MECHANISM FOR ALL DENSITIES 


Test No. 

Void Ratio 
e 

ave. 

Water 
Content, w 
(%) 

Predicted(q^) 
Measured (q^) 

D/B=5 

D/B=10 

D/B=15 

A-l* 

1.030 

1.66 

3.12 

2.95 

2.85 

B-l* 

0.745 

1.60 

1.25 

1.43 

1.30 

B-3 

0.684 

1.75 

1.15 

1.20 

1.25 

B-2 

0.570 

1.40 

1,03 

1.02 

1.02 


*Average of (2) tests. 

soils. However, the general validity of this semi-empirical procedure has not 
yet been established. 

SUMMARY 

Comparisons of unit penetration resistance values measured in controlled 
laboratory tests and predicted by the proposed analytical method have been 
made for both cohesionless (Monterey Sand No. 0) and cohesion-friction (Lunar 
Soil Simulant No. 2) soils. The agreement between measured and predicted 
values for the sand and dense simulant were quite good, suggesting that the 
proposed analytical method may be used confidently for the prediction of the 
static penetration resistance of relatively incompressible soils. A method 
has been suggested for determination of bearing capacity factors N* and N*^ 
for compressible soils. However, the general validity of this method has 
not yet been established. 
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Unit penetration resistance, q f - psi 



FIG. 6.19 COMPARISON OF MEASURED AND PREDICTED 
RESISTANCE CURVES FOR DENSE LSS No. 2 
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FIG. 6.20 FAILURE MECHANISM FOR 
COMPRESSIBLE SOILS 

{After Vesic, 1963) 
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CHAPTER SEVEN 


DETERMINATION OF IN-SITU SOIL PROPERTIES 


INTRODUCTION 

Because the static penetration resistance to a given penetrometer- is- con- 
trolled by the soil shear strength and compressibility characteristics, it 
should be possible to deduce the in-situ properties from the results of static 
penetration tests. A procedure for doing this is described in this chapter. 
This procedure is illustrated by application to laboratory model tests, to 
penetration data obtained from the literature, and to penetration data for 
the lunar surface obtained by the Apollo 15 self-recording penetrometer and 
the Soviet Lunar Rover Lunokhod-1 . 


PROCEDURE 

As discussed in Chapter Four, the ultimate, penetration resistance (q^) 
may be calculated using the following equation: 


q. = cN £ + Y BN £ 

H f c 4 c 's Yq Y<1 


(7.1) 


where c = cohesion 

y = unit weight of soil (pg) 
s 

B = width or diameter of penetrometer base 
£ c , £^ q = shape factors 

N c , N q = bearing capacity or ultimate resistance factors, f^, 
(4>, 6/<j>, a, D/B) 

<{> = angle of internal friction (triaxial residual values) 
6 = friction angle between penetrometer base and soil 
a = base semiapex angle 

D/B = relative depth of penetrometer base , 

Cohesionless Soils (c=0) 

For cohesionless soils (c=0) , Equation 7.1 simplifies to: 


or, in dimensionless form: 


q = Y BN £ 

4 f 's yq Yq 


q f /(Y s B) « N Yq £ Tq 


(7.2) 


(7,3) 
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For a given static penetration test, the following information will be known: 
Cone resistance versus penetration depth or, in dimensionless form. 


( 1 ) 

( 2 ) 


the resistance parameter q / (y B) versus relative depth (D/B) . 

1 S 


Cone semiapex angle (a) and cone roughness (6/4>) . (Cone roughness 
can be approximately estimated from previous experience, from direct 
shear interface test results, and from the surface characteristics 
of the penetrometer material. 

It can be shown that strength properties calculated from Equation 7,1 are 

not very sensitive to soil unit weight, hence any reasonable value can be 

assumed for the calculations. Therefore, the value of N £ can be calculated 

yq yq 

from Equation 7.3. In functional form, N £yq ma y t> e expressed as follows: 


N yq c yq = n * 9 6/c|) ’ a ’ D/B) 


(7.4) 


The only unknown in Equation 7.4 is the angle of internal friction (cj>) , and 
thus (f) may be estimated from this expression for known values of a, 6/<j), and 
D/B using the following procedure: 

(1) Assume a value of (p 

(2) Enter figures in Appendix A relating to cp for the given para- 
meters (a, 6/(j), D/B) and find corresponding value of 

(3) From Equation 4.32 determine the value of g for the assumed (p and 
known relative depth 

(4) Repeat steps 1 through 3 for other values of friction angle (cp) 

(5) Prepare a plot of N £ versus (p (for the specific values of D/B) 

i H I H 

(6) Having already calculated the actual value of Nyq^yq from Equation 7.3, 
enter the plot prepared in step 5 and determine the correct value of (J> 

This procedure may be repeated for different values of D/B to obtain the varia- 
tion of <j) with depth, 
is the assumption that the at rest lateral earth pressure coefficient, K = l-sin<f>. 
For other values of K new values of N should be determined using Equation 4.9. 

yq 

Cohesive Soils (j>=0) 

For purely cohesive soils (4>=0) , Equation 7.1 simplifies to: 


Implicit in the use of the values of N from Appendix A 


q, ■ cN 5 + y D 

H f c c 's 


or, in dimensionless form: 


q £ /c = N ? + 

H f c c 


VI 


(7.5) 


(7.6) 
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As mentioned previously, N £ may be expressed in functional form as: 

N c?c = F(4>, 6/<f>, a, D/B) (7.7) 

For (j>=0, Equation 7.7 simplifies to: 

N c 5 c = F*(a, f c , D/B) (7.8) 

where f is equal to c /c (ratio of adhesion to cohesion), 
c a 

Utilizing depth factors, Equation 7.8 may be rewritten as: 


N 5 = d *N *(a, f )? 

c c c c c c 


(7.9) 


where N* (a, f c ) = bearing capacity cohesion factor for shallow foundations 

obtained from Fig. 2.12a 

£ c = shape factor for cohesion term (equal to 1.2 for <f)=0) 
d c = depth factor for cohesion term 

= loO + b/ ' d'* ' + ' 5 ' o"" 6 (f rom Equation 4.19a) 

Thus Equation 7.5 may be rewritten as: 


q c » cd N *E, + y D (7.10) 

^f c c c s 

From static penetration data, q^ is known for the range of relative 
depth tested. Consequently, from Equation 7.10, cohesion may be calculated 
directly for known values of q^, D/B, a, and f c . This procedure may be 
repeated for other values of D/B to obtain the variation of cohesion with 
depth. 



For cohesion-friction soils, the problem of determining in-situ soil 
properties from the results of static cone penetration test is more com- 
plicated for the following reasons : 

(1) There are two unknowns, c and <j), instead of one 

(2) Most c-cj> soils are more compressible than most cohesionless soils 

For c— 4> soils, the ultimate static cone penetration resistance may be written 
in dimensionless form as follows: 

q £ /(Y s B) - (c/Y s B)N c 5 c + N Yq ? Y<i (7.11) 
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For static penetration tests performed with A given cone there are many 
combinations of c and <f> which satisfy Equation 7.11 for a given value of 
q^/ (Y g B ) • The procedure for developing the c-<j> relationship for given values 
of a, S/4>, and D/B is as follows: 

(1) Assume a <p value 

(2) Enter figures in Appendix A relating and .N to <j> for the 

given parameters (a, 8/<p, D/B) and find the corresponding values 

of N and N 
c yq 

(3) From Equations 4.24a and 4.32 determine the values of £ c and 

£ for the assumed d> and known D/B 
Y q 

(4) Substitute the values obtained in steps 2 and 3 into Equation 7.11 
and solve for cohesion 

(5) Repeat steps 1 through 4 for other values of <f> 

(6) Plot the resulting c-<p relationship 

Steps 1 through 6 may be repeated for other values of relative depth to 


obtain the variation of the c-<p relationship with depth. 

It is obvious that in order to estimate c and <j) separately, more 

-t 

information must be provided. The most convenient way to obtain this in- 


formation is to perform an additional penetration test in the same soil 


deposit using a cone of a different size but having the same values of a. and 
6/<f> as the first. With the penetration resistance data from two sizes of 
cones at the same relative depth (D/B) , the following relationships may 
be written from Equation 7.11: 



1 


(q-/Y B). - N £ 
t 'g/l Yq Yq 

N £ 
cc 


(7.12) 



2 


( w>2 - h yAc 

N c«c 


(7.13) 


where the subscripts 1 and 2 indicate data from the small and large cones 
respectively (B2=kB^) . It is of interest to note in the above equations 
(7.12 and 7.13) that the unit penetration resistance for a large cone will 
be less than for a small cone at a given depth D, Assuming that there is 
no significant variation in cohesion over the depth range of D=mB^ to kmB^ 
(m=D/B) the following relationship may be obtained from Equations 7.12 and 7.13: 
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(7.14) 


B_ (q./y B). - N t 
2 _ H f 's 1 yq yq 


B 


1 (q./y B) - N £ 

H f 's 2 Yq yq 


An expression for ^yq^yq ma ^ Be °b ta i ne d fr° m Equation 7.14 as follows: 


Vt, ‘ 


B 2^ q f //y s B ^2 ~. B l^ q f^ Y s B ^l 
B 2' B 1 


(7.15) 


The value of ^yqSyq ma y calculated directly from Equation 7.15. The 
procedure for the determination of d) from a known value of N £ was 

yq yq 

given previously for the case of cohesionless soils. From Equations 7.12 
and 7.13 the following relationship may be obtained for the determination 
of cohesion: 


Y s A V Y » B) 

N £ A(l/B) 

c c 


(7.16) 


or, in dimensionless form: 


Y s B l N £ 

c c 


(7.17) 


where k = ratio of B^ to B^ 

A = difference for given D/B 

It should be noted that the value of cohesion obtained by two-size cone 
method using Equations 7.16 and 7.17 for a given relative depth represents 
an average value of cohesion for depths in the range of H-mB^ to kmB^. 

It is also possible, if the penetration data from the two-size cone 
method is presented in the form of q^ versus D, to estimate the soil friction 
angle and cohesion at a giveii depth. The following relationships may be 
written from Equation 7.11: 


fcj (q f / T 3 B)i - 

M 1 " 'Wi 


(7.18) 
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2 



(7.19) 


where the subscripts 1 and 2 indicate data from the small and large cones 
respectively. The following relationship may be formulated from Equations 
7.18 and 7.19 by noting that £ c is independent of relative depth: 


f2 . fW>l - ^VcVl l (tl c>2 

h [Cq £ /Y S B > 2 - (N Yq 5 Yq )J (Vl 


(7.20) 


In order to determine <J> from Equation 7.20, the following procedure may 
be used: 

(1) Assume a value of <p. 

(2) Determine values of N c and N corresponding to the assumed value 
of <{> from the figures in Appendix A for the given paramenters 

(a, D/B). It must be remembered that the values of (D/B^) 

and (D/B^) are not the same. 

(3) Determine the values of £ from Equation 4.32 corresponding to 
the assumed value of <p and known D/B. 

(4) Calculate the ratio B 2 /B^. 

(5) Repeat steps 1 through 4 for other assumed values of <j>, 

(6) Prepare a plot of versus $ • 

(7) Enter plot with known value of B 2 ^i ani * the correct value 

of cf>. 

This procedure may be repeated for different values of D to determine the 
variation of <j) with depth. It should be noted that because the contribution 
of cohesion relative to friction^-surcharge decreases with increasing depth, 
the method is not expected to be a sensitive technique for determination of 
cohesion at large depths. 

It should be noted that same method could be applied even if there are 
only data available for one size cone. In this case penetration resistance 
values (q^) corresponding to two different penetration depths should be 
considered. The ratio (B 2 /B^) in Equation 7.20 will be equal to unity. 
Application of the method for this special case is only valid for homogeneous 
soil deposits. 
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From Equations 7.18 and 7.19, the following relationship may be obtained 
for the determination of cohesion: 


c = 


A(q f ) - T s A(BN Yq C Yq ) 

W V 


(7.21) 


The values of N , N , and £ obtained in order to calculate the differences 

c yq yq 

(A) in Equation 7.21 must, of course be selected at the appropriate relative 
depths (D/B) , 

It is also possible to estimate c and cj) separately at the same penetra- 
tion depth (D) if the two cones used are not only different in size but also 
have different values of a and 6/(|>. The procedure for determining $ is the 
same as described previously using Equation 7.20. However, in step 2 of the 
procedure, the correct values of a and 6/<f> must be used. Cohesion may also be 
determined using Equation 7.21. 

The procedures for determining the friction angle and cohesion of relatively 
homogeneous soils using penetration data obtained by the two-size cone method 
needs experimental verification. Also, the bearing capacity factors employed 
in these procedures are valid only for general shear failure conditions. For 
compressible soils the bearing capacity factors will be smaller than those 
corresponding to relatively incompressible soils. Consequently, the use of 
general shear bearing capacity factors in problems involving compressible soils 
will lead to an underestimation of shear strength parameters. 


APPLICATIONS 

Introduction 

The procedure for cohesionless soils is illustrated by application to 
laboratory penetration tests data for Monterey Sand No. 0 and to data presented 
by Green (1970) for Yuma and Mortar sands. The procedure for cohesion-friction 
soils has been applied to laboratory penetration data for LSS No. 2, to lunar 
surface penetration data obtained during the Apollo 15 mission, and by the 
Soviet Lunar Rover Lunokhod-1 . Where possible, comparisons are made between 
actual (measured) and deduced in-situ soil properties. 

Cohesionless Soils 

Monterey Sand No. 0 

The procedure for the determination of the in-situ properties of 
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cohesionless soils was applied to the results (see Chapter Six) of static 
wedge and cone penetration tests in Monterey Sand No. 0, The results for 
wedges are summarized in Table 7.1 and presented in Fig. 7.1. The results 
for cones are summarized in Table 7.2 and presented in Fig. 7.2. As can be 
seen in Figs. 7.1 and 7.2, the agreement between the predicted and actual 
values of soil friction angle is quite good. 

Yuma and Mortar Sands 

The procedure was applied to measured average penetration resistance 
values presented (see Fig, 7.3) by Green (1970) for Yuma and Mortar sands 
using rough flat ended penetrometers (plates) in the range of 5.08 cm to 
60.69 cm in diameter. The gradation curves for Yuma and Mortar sands are 
given in Fig. 7.4. Results are summarized graphically in Fig. 7.5. It can 
be seen that predicted and measured values of <j) agree well. 

Cohesion-Friction Soils 

Lunar Soil Simulant No. 2 

The procedure for the determination of the in-situ properties (cohesion 
and friction angle) of cohesion-friction soils was applied to the results of 
static cone penetration tests in dense Lunar Soil Simulant No. 2 (see Fig. 
6.19). Because data were available for one-size cone only penetration 
resistance values (q^) corresponding to relative depths of 5, 10, and 15 
were considered. This was justified because it was known that the soil pro- 
file was homogeneous. By the application of Equations 7.20 and 7.21 and 
the previously described procedure, strength parameters of c=0.137 psi. 

(0.945 kN/m 2 ) and <j>= 46.6° were obtained. Directly measured values of strength 
parameters for the simulant, corresponding to e=0.57 and w=1.4% were 
c=0.121 psi. (0.834 kN/m 2 ) and <j>=47.0° (triaxial) . It can be seen that 
agreement between predicted and measured values of both c and <}> are excellent. 

Apollo 15 Soil Mechanics Experiment 

The Apollo 15 mission to Hadley Rille provided for the first time in the 
U. S. Lunar Program, quantitative measurement of forces of interaction 
between a penetrometer and the lunar soil. Quantitative analysis of the in- 
situ mechanical properties (c and <J>) of the lunar soil can be made using the 
data obtained by the self-recording penetrometer (SRP) shown in Fig. 7.6. 

The SRP was used tp obtain penetration versus force data in the upper part 
of the lunar surface. The Apollo 15 SRP could penetrate to a maximum depth 
of 76.0 cm (30.0 in.) and could measure penetration force to a maximum of 111 N 
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Table 7.1 SUMMARY OF PREDICTED FRICTION ANGLES 
OF MONTEREY SAND NO. 0 


Depth Range 

Void 

Ratio, e 

Predicted 4> -deg. 

(1) 

. 

(2) 

(3) 

(4) 

(5) 

D/B<5 

0.728 

39.5 

39.5 

39.2 

39.7 

39.5 

O'. 665 

41.2 

41.7 

41.2 

42.0 

42.1 

0.599 

44.5 

44.3 

43.8 

44.2 

45.3 

D/B<10 

0.728 

40.0 

40.1 

39.5 

39.7 

40.1 

0.665 

41.7 

42.7 

41.5 

42.0 

43.7 

0.599 

44.5 

45.5 

44.9 

45.1 

46.1 

D/B<15 

i 

0.728 

- 

40.2 

- 

39.7 

40.1 

0.665 

42, 2 

43.7 

42.0 

42.1 : 

44.0 

0.599 

44.8 

46.5 

45.2 

-i 

44.9 

46.1 


(1) a - 90°, S/4>= 0.5; (2) a = 15% 6/<f>= 0.5; 

(3) a = 15% 6/<f>= 0.9; (4) a = 30% 6 / 4 >= 0.9; 

(5) a = 30% 6/<j>= 0.5 (B-0.8 in. Wedges) 


Void Ratio 
e 

Predicted* 
<j)-deg . 

Measured** 

^peak - de §- 


0,728 

39.8+0.3 

39.7 

39.2 

0.665 

42.2+0.6 

42.6 

41.8 

0.599 

45.0+0.5 

45.4+0.2 

44.5+0.2 


* For all depths and penetrometer configurations 
**Triaxial friction angles 
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FIG. 7.1 COMPARISON OF MEASURED AND 
PREDICTED FRICTION ANGLES 
FOR MONTEREY SAND No. 0 
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Table 7.2 SUMMARY OF PREDICTED FRICTION ANGLES 
OF MONTEREY SAND NO. 0 


Depth Range 

Void 
Ratio , e 

Predicted cj> -deg. 

(i) 

(2) 

(3) 

(4) 

(5) 

D/B<5 

0.728 

39.1 

41.7 

39.5 



0.665 

40.4 

- 

42.0 

- 

- 

0.659 

& 

42.6 

- 

- 

- 

0.599 

43.8 


44.9 

44.8 

44.8 

0.575 

- 

45.0 

- 

- 


d/b<lo 

0.728 

39.1 

40.2 

39.1 

- 


0.665 

41.0 

- 

42.0 

- 


0.659 


42.8 

- 

- 

- 

0.599 

44.2 

- 

44.9 

45.2 

45.4 

0.575 

- 

46.3 

- 

- 


D/B<15 

0.728 

39.1 

40.0 

? 

- 

- 

0.665 

41.7 

- 

? 

- 

- 

0.659 

- 

43.2 

- 

- 

- 

0.599 

44.5 

- 

45.1 

45.3 

45.7 

0.575 

- 

45.5 

- 


- 


(1) a = 90°, 6/<j> - 0.5; (2) a = 15°, 6/<j> = 0.5 

(3) a = 15°, 6/4> = 0.9; (4) a = 15°, 6/4 = 0.3 

(5) a = 30°, 6/4 = 0.3 (B = 0.8 in. Cones) 


Void Ratio 
e 

Predicted* 
<j>-deg . 

Measured** 



0.728 

39.7 

39.7 

39.2 

0.665 

41.5 

42.6 

41.8 

0.659 

42.9 

42.9 

42.0 

0.599 

44.9 

45.6 

44.7 

0.575 

45.6 

46.8 

45.6 


* For all depths and penetrometer configurations 

**Triaxial friction angles 
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Measured , <f> - degrees 


FIG. 7.2 COMPARISON OF MEASURED AND 
PREDICTED FRICTION ANGLES 
FOR MONTEREY SAND No. 0 
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Relative Depth D /g 

MORTAR SAND 

FIG. 7.3 PENETRATION CURVES FOR YUMA AND MORTAR SANDS 

(After Green, 1970) 
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FIG. 7.4 GRADATION AND CLASSIFICATION DATA FOR YUMA 
AND MORTAR SANDS (After Green, 1970 ) 
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FIGt 7.5 COMPARISON OF MEASURED AND PREDICTED FRICTION ANGLES FOR YUMA AND MORTAR SANDS 

(Data after Green, 1970) 









(25 pounds). The record of each penetration was scribed on a recording drum 
contained in the upper housing assembly. Detailed information about the SRP 
is given by Mitchell et al (1972a) . 

As reported by Mitchell et al, (1972a), because of the tendency of the 
lunar reference plane to ride up on the penetrometer shaft precise values 
of penetration were not known for the penetration tests on Apollo 15, and 
the exact shape of the force versus depth curve was not obtained. None- 
theless, estimates of the penetration were possible as indicated in Table 7.3. 

Table 7.3 SUMMARY OF APOLLO 15 CONE PENETRATION TEST 
RESULTS FOR 0.5 SQ. IN. CONE 
(after Mitchell et al, 1972a) 


Location Near 
Traverse Station 8 

Penetration 
q^=50.0 psi. , 

at 

(cm) 

Relative 
Depth (D/B) 

Adjacent to soil 




mechanics trench 

8.25 


4.06 

In LRV track 

5.25 


2.58 


The c-<p relationship for the lunar surface material may be deduced by 
using the procedure described earlier and the information listed below: 

(1) Roughness (6/4)) is approximately equal to 0.5. (Based on the 
results of friction measurements between a ground basalt lunar 
soil simulant and hard-anodized aluminum similar in roughness to 
that used for the SRP cones) . 

(2) Unit weight of soil y -1.8 g/cm 3 (under earth gravity). 

s 

(3) Ultimate resistance q^=50.0 psi. (34.5 N/cm 2 ). 

(4) Measured values of relative depth. (Table 7.3 indicates that the 
relative depths (D/B) for these penetration tests fall in the 
range of about 2.5 to 4.1). Thus the value of D/B=3 may represent 
the actual conditions reasonably well. 

The results of these calculations are shown in Fig. 7.7. Another c-(j) 
relationship was established during the Apollo 15 mission by excavating a 
trench and failing the wall of the trench" by ar krrown boundary loading. The 
c— 4> relationship for the soil mechanics trench has also been plotted on 
Fig. 7.7. The intersection of this curve and the c-(}> curve for the penetra- 
tion tests provide values of c and tj) which simultaneously satisfy the results 
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FIG. 7.7 PROPERTIES OF LUNAR SURFACE DEDUCED . 

FROM APOLLO 15 PENETRATION TEST RESULTS 
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of both the trench and penetration testa. For D/B=3, the required cohesion 
is 1.1 kN/m 2 (0.16 psi.) and the angle of internal friction is 46.0°. It 
should be remembered that this value is expected to be close to the residual 
value of friction angle. The value obtained agrees closely with that 
obtained by comparison of the observed penetration behavior with that of 
terrestrial simulants. 

Detailed analysis of Apollo penetration tests are presented in the 
Preliminary Science Reports for Apollo 15 and Apollo 16 (Mitchell et al, 

1972a and Mitchell et al 1973). 

Lunokhod-1 

The procedure for deduction of cohesion and friction angle from the 
results of static cone penetration tests was applied to penetration resistance 
data for the lunar surface obtained by the Soviet Lunar Rover Lunokhod-1 which 
were obtained in the western side of Mare Imbrium. Data were obtained from 
the report of Leonovich et al (1971) . 

Investigations with the Lunokhod-1 were made using a special instrument 
by penetrating and rotating a conical'-bladed punch in the ground. With the 
penetration of the punch it was possible to determine ground carrying 
capacity*, and with rotation of the vanes the resistance to rotational shear 
(torque resistance**) was determined. During the three months of Lunokhod 
operation there were 327 measurements of mechanical properties by means of 
penetration of the conlcal-bladed punch at different ground locations; e.g. , 
craters, concentrations of rocks, horizontal and sloping surface sectors. 

Summaries of the data provided by Leonovich et al (1971) are given in 
Tables 7.4 and 7.5. Fig. 7.8 shows a histogram and curve of the statistical 
distribution of ground carrying capacity and torque resistance along one of 
the trajectory segmants. 

Analysis of vane shear tests 

Usually interpretation of vane shear test results is based upon the 
assumption of a uniform distribution of shear resistance on the whole surface 
of the rotating cylinder, including the two ends. The same assumption appears 
to have been made by Leonovich et al (1971) leading to an estimation of the 


* Ground Carrying Capacity was defined as the specific resistance on the punch 
for a depth of ground penetration equal to the height of the punch. 

**Torque Resistance was defined as required cohesion neglecting the contri- 
bution of friction to vane resistance. 
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Table 7.4 SUMMARY OF TECHNICAL DATA FOR CONE-VANE 
PUNCH (after Leonovich et al 1971) 


Maximum force on the penetrometer 


196 N 

Depth of penetration 


50-100 mm 

Angle of punch rotation 

up 

to 90 deg. 

Maximum torque on punch 


0,5 kgm 

Diameter of cone base 


50.0 mm 

Diameter of blades 


70.0 mm 

Height of cone-blade punch 


44.0 mm 

Blade thickness (square cut edge)* 

1.1 

to 1.6 mm 

Surface roughness of cone* 

0.3 

to 0.5 

*Estimated from the photographs 




Table 7.5 SUMMARY OF DATA FOR THE MECHANICAL 
PROPERTIES OF LUNAR GROUND 
(after Leonovich et al, 1971) 


Range of Values 


Carrying capacity 

19.6 to 98.0 kN/m 2 

Torque resistance 

1.96 to 8.82 kN/m 2 

Density of soil 

1.5 to 1.7 g/cm 3 

Highest frequency values 


Carrying capacity 

33.3 kN/m 2 

Torque resistance 

39.6 kN/m 2 
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Carrying capacity 


(a) 


1.31 5.40 9.69 kN/m Z 



Torque Resistance 


(b) 

FIG. 7.8 HISTOGRAM AND CURVE REPRESENTING DISTRIBUTION 
OF GROUND CARRYING CAPACITY AND GROUND RESIS- 
TANCE TO TORQUE (After Leonovich et al, 1971) 
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highest frequency cohesion (or torque resistance) value of 39.6 kN/m 2 (0.683 
psi) . The estimated range of cohesion from Surveyor and Apollo data was only 
0.345 to 1.38 kN/m 2 . The high cohesion value estimated from the Lunokhod-1 
data may be due to the fact that the contribution of friction to torque 
resistance was ignored by Leonovich et al, (1971). Even if it was assumed 
that a cylindrical failure pattern was actually obtained, the distribution 
of stresses around the circumference could not be uniform for soils in which 
there is an appreciable contribution of friction to the resistance, A uniform 
shear strength distribution is correct only for 4>=0 soils, where shear strength 
is independent of normal stress. 

Farrent (1960) derived the following expression for torque in c~4> soils 
by assuming that shear strength, s q , develops at both ends and immediately 
behind the blades: 


m 4 _ 3 . nr 2 L 

T = -r irr s + r- — t“ 
3 o 2tan(p 


s (e 


4rr 

n 


tan<p 


- 1 ) 


(7.22) 


where T = torque 

r = radius of vane 

n = number of blades 

L = length of vane 

s q = shear strength of soil 
<f> = angle of internal friction of soil 

When the contribution of friction is considered, the minumum value of 
shear strength (s q ) can be calculated by the following relationship (after 
Cox, 1967): 


s q = c(l + sin<J>) 


(7.23) 


For the analysis of Lunokhod data (n=4 and L=1.25r) the following 
relationship can be written: 


T - irr’cU+sin*) + -i) 


or, in dimensionless form: 


r 3 c 


4 . 5(e 1,tm> -1) 

3 2tantp 


(l+sin<}>) 


(7.24) 


(7.25) 
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It can be shown that for cf>=0 soils Equation 7.25 reduces to: 


T 

r 3 'c 


= fir 


+ fir 


23 tt ~ 
6 " 


12.0 


(7.26) 


The same equation can be obtained by considering the development of uniform 
shear strength (cohesion) on the both ends of the vane and on the cylindrical 
failure surface. 

Values of the dimensionless parameter T/ (r 3 c) for different values of <f> 
were calculated from Equation 7.25 and are summarized in Table 7.6. 


Table 7.6 VALUES OF DIMENSIONLESS PARAMETER T/(r 3 c) 
FOR DIFFERENT VALUES OF <f> 


<f>(deg.) 

20 

25 

30 

35 

40 

45 

50 

a 

CO 

u 

h“ 

25.4 

31.3 

39.6 

52.4 

70,5 

102 

161 


It can be seen that T/ (r 3 c) increases with increasing 4> for c— <J> soils indicat- 
ing that the necessary cohesion should decrease with increasing <j) for a given 
T (torque) value. Consequently , the ratios of T/ (r 3 c) corresponding to e-<|> 
and cf>=0 analysis can be determined from the following relationships: 


lr 3 cj (<t>=0) 


23tt 

6 


r 3 c (c-<J>) 


■ 5(e ir CaP <(, _ x) 

3 2tan<|) 


(l+sin())) 


or 


(c) 

(c) 


(d?=o) = 

(c-4>) 


4rr 5(e 7rtan<t> -1) 
.3 2 tantj) 


(l+sin<|)) 


23tt 

6 




(7.27a) 


(7.27b) 


The calculated values are summarized in Table 7.7 below: 
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Table 7.7 THE RATIOS OF COHESION VALUES CALCULATED 
FROM <f>=0 AND c-<f) ANALYSIS 


<)>(deg.) 

20 

25 

30 

35 

40 

45 

50 

C (d)=0) 

C (c-cj>) 

2.11 

2.61 

3.29 

4.36 

5.85 

8.46 

13.4 


The value of cohesion c=39.6 kN/m 2 (0.683 psi) reported by Leonovich et al 
(1971) corresponds to a (p=0 analysis. Corresponding c-<|> relationships were 
calculated using Table 7,7 for c-cf> analysis, and are shown in Fig. 7.9. It 
can be seen that for <f>=45° the calculated value of cohesion from c-<j> analysis 
is 8.5 times smaller than the value calculated from (f>=0 analysis. 

The variation in c— <J> relationships for different ground locations may be 
determined using the range of torque resistance values (1.96 to 8.82 kN/m 2 ) 
reported for these locations. The corresponding c-<f> relationships have been 
calculated and are shown in Fig. 7.10. The c-<f> relationship having the 
highest frequency is also plotted in Fig. 7.10. It may be concluded on the 
basis of these data that the mechanical properties of the lunar ground can vary 
over an extremely wide range. 

Analysis of cone penetration tests 

The cone penetration data from Lunokhod-1 also provide c-cf> relationships 
for the lunar soil. The penetration resistance of the soil to the cone-vane 
punch has two components; (1) resistance due to end bearing of vane blades, and (2) 
resistance due to the cone. The largest contribution to penetration resistance 
comes from the cone, because its bearing area is much larger than the bearing 
area of the vane blades. The highest frequency value of ground carrying 
capacity (q^) is 33.3 kN/m 2 . If it is assumed that the contribution due to 
end bearing of the vane is negligible and the soil disturbance due to the 
penetration of the vane is insignificant, for a=30° , 6/<J>=0.5 (assumed) and 
D/B=0, the calculated c— cj> relationship is as shown in Fig. 7.11. 

The effect of the end resistance of the vanes may be evaluated by 
considering the blades of the vane as two strip footings with a width of 1.5 mm 
and length of 70 mm. For a penetration of 44 mm, the relative depth (D/B) is 
approximately 30. For comparison, the calculated c-<j> relationship consider- 
ing both cone and vane resistance is also plotted in Fig. 7.11. 

The effect of cone roughness must also be considered. Calculations for 
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FIG. 7.9 HIGHEST FREQUENCY C - <f> RELATIONSHIP OBTAINED 
FROM THE ANALYSIS OF VANE SHEAR TESTS 

{Dofo from Leonovich et ol, 1971) 
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Cohesion, 


IO.O, 



i 

FIG. 7.10 RANGE OF c-<£ RELATIONS OBTAINED FROM THE ANALYSIS 
OF VANE SHEAR TESTS (Data from Leonovich et al, 1971) 

>■ i 
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Cohesion 


without considering 
end bearing of vane blades 


'$/<£ * 0.3 (assumed) 


by considering / y \ 

end bearing of vane ■/ \ \ 

blades, ( 8 /$ = 0.5) 

s 0,5 (assumed) 


d/ b = 


Angle of internal friction, </)- degrees 


FIG. 7. II EFFECTS OF CONE ROUGHNESS AND END BEARING OF 
VANE BLADES ON SHEAR STRENGTH PARAMETERS OF 
LUNAR SURFACE (Data from Leonovich etal, 1971} 



Fig. 7.11 were based on the assumption that the cone roughness (S/4 r) was 
0.5. For a smooth cone surface roughness will be close to 0.3, The cal- 
culated c— relationship for 6/4>=»0 . 3 (neglecting end bearing of vane) is also 
shown in Fig. 7.11. It can be seen that the curve for 6/(f>=0.5 with end 
bearing of vanes neglected represents a good average of the three alternatives 
considered. 

The variation in c-cf> relationships for different ground locations may be 
determined using the range of carrying capacity (19.6 to 98.0 kN/m 2 ) reported 
for these locations. The corresponding c-$ relationships have ijeen calculated 
and are shown in Fig. 7.12. The e-<j> relationship having the highest fre- 
quency (q f = 33.3 kN/m 2 ) is also plotted in Fig. 7.12. 

Comparison of c— relationships 

A comparison of the c-4> relationships obtained from two different analyses 
(vane shear and cone penetration) is shown in Fig. 7.13. Considering the 
different failure mechanisms involved, the c-<p relationships are quite 
comparable. The test results for the more dense soils (higher values of tj)) 
show higher values of strength parameters for the vane test results than for 
the penetration tests. This may be because the vane tests were done at 
somewhat greater depths, and also becuase the analysis of vane data may be 
expected to yield values of <t> close to the peak; whereas, the cone results 
reflect more closely the residual strength parameters. The difference between 
peak and residual values for most soils generally increases with increasing 
density (and therefore increasing <f>). Further, the c-<$> relationships developed 
from the cone-vane penetration data reflect the assumptions made regarding 
appropriate value of 6/4> for the cone (a value of 0.5 was used for the curves 
shown in Fig. 7.13) and the contribution of the vane (neglected for the curves 
shown in Fig. 7.13). 

Change in shear strength properties with depth 

Leonovich et al (1971) present separate penetration data in the form 
of force versus penetration for different ground sectors as shown in Fig. 7.14. 
Separate analysis of each individual curve can be made, and the change in 
shear strength properties with depth can be studied . This can be done only 
if the change in bearing capacity factors with depth are formulated. As 
discussed in detail in Chapter Four, bearing capacity factors can be cal- 
culated from Equations 4.6 and 4.9 as functions of relative depth (D/B) . 

As an example, the penetration data for horizontal ground sector (curve 1 
in Fig. 7.14) were analyzed. The penetration data obtained from Fig. 7.14 
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Cohesion 



FIG. 7.12 RANGE OF SHEAR STRENGTH PARAMETERS OBTAINED FOR 
LUNAR SURFACE FROM CONE PENETRATION TEST RESULTS 

(Data from Leonovich et al , 1971) 
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FIG, 7.13 COMPARISON OF fc-<£ RELATIONS OBTAINED 
FROM CONE PENETRATION AND VANE SHEAR 
TESTS 
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Load on the punch, Q- Newtons 
0 40.2 80.4 120.6 160.8 201.0 



FIG. 7.14 PENETRATION DATA OBTAINED BY LUNOKHOD-1 
FOR FOUR CONDITIONS IN THE WESTERN PART 
OF MARE IMBRIUM 

(Data from Leonovich, et al. 1971) 
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are summarized in Table 7.8 for different relative depths. 


Table 7.8 PENETRATION DATA FOR LUNAR SURFACE AT A 
HORIZONTAL GROUND SECTOR 
(after Leonovich et al, 1971) 


Relative Depth 
(D/B) 

Force 
Q (kg) 

Unit Penetration 
Resistance, 

kg/ cm 2 

kN/m 2 

0.00 

7.2 

0.367 

36.0 

0.10 

8.0 

0.408 

40.0 

0.25 

9.8 

0.500 

49.0 

0.50 

12.8 

0,653 

64.0 

0.70 

15.8 

. i 

0.806 

79.0 


For a given relative depth (D/B) using the procedure outlined in this 
chapter, corresponding c— 4> relationships were calculated. The calculated 
relationships are shown in Fig. 7.15. It can be seen that there is a slight 
increase in shear strength parameters with depth. This may be an indication 
of an increase in soil density with depth. 

SUMMARY 

Procedures for the determination of in-situ soil strength parameters 
from the results of static penetration tests in cohesionless, purely cohesive, 
and cohesion-friction soils have been presented. These procedures have also 
been illustrated by application to laboratory test results, to penetration 
data obtained from the literature, and to penetration data for the lunar 
surface obtained by the Apollo 15 self-recording penetrometer and the Soviet 
Lunar Rover Lunokhod-1. Comparisons of measured and predicted strength para- 
meters have been made where possible. These comparisons indicate very good 
agreement between measured and predicted values. 
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FIG. 7.15 CHANGE IN SHEAR STRENGTH PROPERTIES 
WITH DEPTH, LEVEL INTERCRATER REGION 
(Dots from Leonovich, et ol 1971) 
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CHAPTER EIGHT 


SUMMARY AND CONCLUSIONS 


GENERAL 

This study was undertaken to develop practical analytical procedures 
which would provide a theoretical basis for the understanding and utilization 
of the results of static penetration tests for the determination of in-situ 
lunar soil properties in support of the Lunar Soil Mechanics Experiment * 

Model tests were used to establish a rational basis for theoretical develop- 
ments and to improve the understanding of penetrometer-soil interaction. 

The effects of penetrometer to soil friction, soil friction angle, base apex 
angle, and relative depth (ratio of penetrometer base depth to the base 
width) on the failure mechanism were investigated. A. failure mechanism has 
been proposed based on the interpretation of features observed in the models 
and the results of previous studies. Knowledge of the failure mechanism 
has permitted the development of new relationships for determining the ultimate 
base resistance which account explicitly fpr such important parameters as base 
apex angle, base roughness, soil friction angle, and relative depth. 

Values of input variables (soil friction angle and penetrometer to 
soil friction angle) for a fine sand (Monterey Sand No. 0) and a silty fine 
sand (LSS No. 2) were determined in order to experimentally evaluate the 
proposed relationships. Carefully controlled laboratory penetration tests 
were conducted on both of these soils. The measured values of input variables 
were used to predict the penetration resistance of these soils in accordance 
with the theory presented » Comparisons were made of the measured and pre- 
dicted penetration resistance values, and agreement was found to be very 
good. 

Analytical procedures for the determination of in-situ soil strength 
parameters from the results of static penetration tests in cbhes ionless , 
purely cohesive and cohesion-friction soils have been presented. These 
procedures have been illustrated by application to (1) laboratory test results, 
(2) penetration resistance data obtained from the literature, and (3) penetra- 
tion resistance data for the lunar surface obtained by the Apollo 15 self- 
recording penetrometer and the Soviet Lunar Rover Lunokhod-1. Comparisons 
of measured and predicted strength parameters have indicated very good 
agreement * 
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FAILURE MECHANISM 

From the results of model tests it was found that a failure surface as 
shown in Fig. 3.11 represents closeiy the actual failure surface associated 
with wedge penetration. A plane shear zone exists adjacent to the base of 
the penetrometer. The dimensions of this plane shear zone depend on both 
the penetrometer to soil friction angle (6) and soil friction angle (<j>) , as 
described by Equation 3.9. 

A logarithmic spiral bounds a radial shear zone to a point of vertical 
tangency, above which the failure surface rises vertically to the ground 
surface for large depths of penetration. For shallow penetration depths, the 
logarithmic spiral breaks out at ground surface before vertical tangency is 
reached. This failure mechanism associated with the penetration of wedges is 
consistent with those suggested also by Biarez et al (1961) and Hu (1965) . 

THEORY 

Equations 4.3, 4.6 and 4.9 describe the resistance to penetration of 
wedge shaped penetrometers in terms of soil friction angle, cohesion, and 
density and base apex angle (2a), base roughness (6/<j>) , and relative depth 
of penetrometer base (D/B) . This solution to the static penetration problem 
provides a generality not previously available. 

For determination of penetration resistance of cones, bearing capacity 
factors corresponding to strip loading must be modified by proper shape 
factors. Equations 4.26 and 4.32 can be used for calculation of these shape 
factors. 

It has been shown that the use of bearing capacity factors formulated 
for general shear failure conditions will cause overestimation of the penetra- 
tion resistance of compressible soils, such as loose to medium-dense silty 
fine sands. As a result of soil compressibility, the shear surface is 
restricted to a smaller zone around the penetrometer tip. A procedure has 
been suggested for determination of bearing capacity factors for compressible 
soils. However the general validity of this procedure has not yet been 
established. 

For resistance of wedge and cone shaped penetrometers the following 
conclusions can be drawn from the theory: 

(1) For rough wedges and cones with obtuse base apex angles, the 

bearing capacity factors are nearly equal to those for plane and 
horizontal contact areas. 
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( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


The bearing capacity factors for rough bases increase with decreasing 
values of base semi-apex angle (a) below approximately 15° . 

For perfectly smooth wedges and cones, the bearing capacity factors 
increase with increasing base apex angle. 

The bearing capacity factors for a given roughness should not be 
estimated by linear interpolation between perfectly smooth and 
perfectly rough values. 

Base roughness has little or no influence on the penetration 
resistance encountered by plane and horizontal contact areas in 
cohesive soils. However, at the surface of a cohesionless soil, a 
smooth flat-ended penetrometer will encounter less penetration 
resistance than one having a rough base. 

For relative depths greater than a critical relative depth depen- 
dent on the particular angle of internal friction, base apex angle, 
and base roughness: 

(a) There will be no change in the bearing capacity cohesion 
factor N c with further penetration. 

(b) The increase in bearing capacity friction-surcharge factor 



with further penetration will be proportional to the increase 


in depth. 


DETERMINATION OF STRENGTH PARAMETERS 

From the results of direct shear interface tests, it is concluded that 
the roughness (<5/<J>) has almost an unique value for a given soil and penetro- 
meter. In other words, 6/<p is not sensitive to soil density. It has also 
been shown that 6/4>=0.3 - 0.5 for a smooth cone and 6 / 4>— 0 . 9 for rough surfaced 
cones „ 

The use of friction angles corresponding to plane strain peak values 
will cause overestimation of the penetration resistance of soils. Because 
the ultimate values of shear stresses control the resistance at large defor- 
mations, friction angles corresponding to ultimate (residual) values of shear 
stresses should be used in the analysis of the continuous static penetra- 
tion problem. It is common practice, in shear testing, however, to obtain 
only peak values of the friction angle and to terminate the test before 
accurate ultimate values are obtained. It has been shown that the use of 
peak triaxial friction angles could cause approximately 20 to 30 percent 
over estimation of the penetration resistance for dense sands. Further, the 
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use of peak values of friction angle for less dense sands will yield correct 
values of penetration resistance, because peak and ultimate values of fric- 
tion angles agree closely for such soils. 

Several examples have been presented to show that predicted and measured 
shear strength parameters of both cohesionless and cohesion-friction soils 
agree well. Therefore, it is concluded that the proposed analytical method 
(Equations 7.1 through 7.21) with the aid of the developed theory can 
satisfactorily be used for the determination of in-situ shear strength 
properties of soils which fail in general shear. 

RECOMMENDATIONS 

It is recommended that further efforts be directed to: 

(1) Verification of the proposed analytical method for a greater 
variety of cohesionless and cohesion-friction soils, 

(2) Further investigation of the effect of soil compressibility, and 
development of an analytical method for determination of penetra- 
tion resistance of compressible soils, 

(3) Development of a method for analysis of the penetration resistance 
of layered soils (in which properties differ greatly among layers) . 
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APPENDIX A 


DERIVATION OF FORMULAS FOR 
BEARING CAPACITY FACTORS 


DETERMINATION OF BEARING CAPACITY FACTOR N 

c 

The failure mechanism and associated free body diagram are given in Fig. 

A.l. According to solutions of equilibrium equations for a weightless body, 

the radial shear zone is defined by a logarithmic spiral and any radial 

plane is a failure plane. Therefore OC is a plane where full mobilization 

of shear, strength takes place, consequently stresses cr and T, can be repre- 

b b 

sented by point B on Mohr's envelope. Stresses and are stresses on 
plane OD. Stress O^ can written in terms of as follows: 


°b " a o + cos* [sin( 2 e ^ ) -sin<f>] 

where = c+a^tan(J) 

By substituting Equation A. 2a into Equation A.l: 

(c+a,tan<j>) 

0,-0 + 7 • [sin(2e+*)-sin(f>] 

b o cos<|) 


(A.l) 


(A. 2a) 


(A. 2b) 


By rewriting Equation A. 2b: 

cos* O +[sin(2e+*)-sin*]c 

rt 8 i Q i M i 

b [cos*-tan*[sin(2e+*)-sin*] ] 


(A. 2c) 


From equilibrium of the logarithmic spiral bounded by OBC: 


EM - 0 
o 


(A. 3a) 


r? r 2 0 

vr-°c 2*- + L cr2d9 -° 


(A. 3b) 


From the general equation of a logarithmic spiral (r=r e ® tan( i > | and by substituting 
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c*Cot <f> - CaCot 8 

FIG. A.1 FREE BODY DIAGRAM FOR DETERMINATION 
OF BEARING CAPACITY FACTOR N c 
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into Equation A. 3b and solving for a c , the following relationship can be 
obtained : 


a 

c 


20 tan$ „ , 

e o cr u + . . 

b tan<p 


c , 20 tan<J) . . 
(e o r -1) 


(A. 4) 


Since OB is a slip plane: 


r » c + a tancb 
c c Y 


(A. 5a) 


Substituting Equation A. 4 into Equation A. 5a: 


T c = (c + a b tan<j>)e 2 ®o tan< ^ 


(A. 5b) 


Substituting the value of from Equatibn A. 2c into Equations A. 4 and A, 5b: 


[l+tan 4 >(sin( 25 +<t>)-sin<J>) ]e 2 ^o tan< ^ c+sin<j> e 2 ® 0 tan< ^ 
[cos4>-tan<j>(sin(2£-kj>)-sin<t>) ] 


(A. 6a) 


[ sin (2£+<J) )-sin<J>] c + cos<}) O q ^ tan( ^ 
°c ~ [cos<(>-tan<Ksin(2£+<j>)-sin$).] e ° 

+ c cot<j> (e 2 ®o tan< ^ -,1) 


(A. 6b) 


From Fig. A.l stresses a and T can be written in terms of O and T as follows: 

P P c c 


a 

P 


[ iI^ +sin(2Y_(f,)1 


T 

C 

cos<j> 


c 

tanc{) 


(A. 7a) 


and 


T 

P 


x . cos (2y-(j>) 
c cos<}> 


(A. 7b) 


By substituting Equation A. 6a into Equations A. 7a and A. 7b the following 
expressions can be obtained for and T^: 
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a 

P 


[I4~sin(j> sin(2y-4>)]c e 2 ^o^ an< ^ 
sin<j> costj) 



[ (sin(2£-k|>)-sin(j>) cos<j) c + c q ] • [l+sin$ sin(2y-4>) ]e 20 o tan< ^ 
cos 2 <J>[cos<}>-tan<f> (sin(2£-kf>)-sin<t>) ] 


(A. 8a) 


T 

P 


[ (sin(2£+cf>)-sin<}Oc*cos<3>+a o ] tan<|> cos (2Y-<|>)e 29 o tan< ^ 
~~~~ cos<f>[cos<}>-tan<}> (sin(2£+<j>)-sin<{>)] 


+ cos(2y-t|))e 29 o tan ^ c 

COS<{> 


(A. 8b) 


From equilibrium of wedge 00' A: 


<T * a p + T p tan¥ (A. 9) 

By substituting values of O and T into Equation A. 9 and recalling that N 

P P c 

is the constant corresponding to cohesion(c), the following relationship is 

obtained: 


N = [lj-sin4>sin(2Y-<t>)]e 29 o tan( ^ cos(2y-<t>) ‘tan^e 290 * 311 ^ 
c sin<|> cos <J> cos(J) 

_ 1 + [sin(2g-Hj>)-sin(i>3 [l+sinfr sin(2Y-<<))3e 29 o tan(t> 

tan<t> « cos 2 (j)[ C Q S (j)»tan4)(sin(2C+^)-sin4))] 

+ tan(|> tan^ cos(2y-<ft) fsin(2Y+4>)-sind>3e 29 o tan< ^ 
cos<f>* [cos<J>-tan<j> (sin(2£+<j>)-sin<j>)] 


However, because it is assumed that no shear stress develops on vertical plane 
CD, angle £ can be taken equal to zero, and Equation A. 10 simplifies to: 

N = [l+sin<t> sin(2Y-<l>)3 •e 29 o tan( ^ > _ 1 

c sin<j>*cos<{> tan<J> 

, cos(2Y-<j>) •tan'i' , e 29 o tan ^ 

+ cos(p (A- 1; 
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DETERMINATION OF BEARING CAPACITY FACTOR N 

■ yq - - 

The failure mechanism and associated free body diagram is given in Fig. 

A. 2. From static equilibrium of the body OCEFG, the sum of the moments about 
point 0 must be zero: 

E M = 0 (A. 12a) 

o 


or 


f, *l, + p :> 

b b si 


“si 


P s2* L s2 


+ P 1 * L pl 


+ Pi * L 2 

2 p 


W- *L n 
1 wl 


(A. 12b) 


Note that moments of F forces are zero because they pass through point 0. 
Determination of 

Refering to Fig. A. 2, the area of logarithmic spiral segment OL^L^ is: 


dA - 1/2 r 2 d0 (A. 13a) 

where r “ r^ e ® tan( ^ (A. 13b) 

By substitution, the previous two equations yield the following relation- 
ships : 

dA = 1/2 r? e 20tan< ^ d0 (A. 13c) 


and 

dW = l/2*y *r 2 e 20tan< ^*d0 (A. 13d) 

s 

The moment around point 0 becomes 

dM - dW*L 

w 

where L = 2/3* r n •e 0tan ^ cos (0-3-0) 

w 1 o 

and by substituting Equation A. 14b and Equation A. 13d into Equation A. 14a: 

dM = l/3*y *ri e 3 ® tan< f > cos (0 -0-0) d0 (A. 15) 

s o 


(A. 14a) 
(A. 14b) 


192 








(A. 16) 


thus. 


W l 1 wl-/° dMol/ 3 ’ Y s- r ’ f 


)0tan<j> 


e" ' cos(0 q — 3 — 0) d0 


If the integral in Equation A. 16 is designated as 


Ig =J* cos(0 o -3-9) d0 


then the following relationship is obtained: 


W l L wl = 1/3 


Determination of I 


0 


By partial integration: 


_ 1 3 0tan4> 

i 0 3tan<(> 6 


U « 

cos(8 o -B-9) |° - 3 ^ J° e 36tan * s i„(0 o -B-9)d6 


Similarly the second term may be partially integrated: 


1 3 0tan<{> rn p Q\ 

3 ^ e co3(8 o -B-6) 


1 3 0tan<(> 


9tan 2 <j) 


sin(0 -3- 
o 


0) 


9tan 2 (() 


u 

j° e 30tancf> cos(0o _ 3 _ e) de 


It can be seen that the last term in Equation A. 19b is the same as Igj 
therefore, the equation may be rewritten as: 


(A. 17) 


(A. 18) 


(A. 19a) 


(A. 19b) 


194 



« + stktf J e -{adgr e!0t “* co.(e 6 -e-e) - 


^ * !0tan * ” i ” ( 9 o- b - 9 >; 


(A, 19c) 


Considering the boundary- values, the following expression for Ig is obtained! 

57 — ST /3tan(J)[e 3 ®o ta,a ^cos3 - cos(0 -0)] + 


9 l+9tan 2 <j> 

[ e 3 0 o tanking + 8in (0 o _6)]| 


(A. 20) 


Determination of PjL^ 

From Fig. A. 2 the following relationships can be written: 


OE - r 2 - r ie 


(A. 21a) 


OH m OE sing - r 1 e 9 o tan ^sing 


(A. 21b) 


HE = OE cosg = r^e 9 o tan ^cosg 


(A. 21c) 


and 


P. = l/2y *0H*HE 
1 's 


(A. 22a) 


It is readily apparent that 


P. = 1/2 y rf sing cosg e z ^o tan< ^ 


(A. 22b) 
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and 


(A. 22c) 


= HE/3 = 1/3 r 1 e®o, tan ^cos3 

and thus. 


P n L 1 ■ 1/6 y r? sin3 cos 2 3 e 3 ®o tan< ^ 
1 pi s 


Determination of P 0 L „ 
t P.£ 


From Fig. A. 2: 


EF = OG-OH = D - r^sin3 e®o tan< ^ 

By introducing, m = D/B, Equation A. 24a may be written as 


EF ■ mB - r^sin3 e®o tan< ^ 


The force is given by: 


P_ = Y * EF*HE 

4 S 

By substituting Equations A. 21c and A. 24b into Equation A. 24c: 


V r i cose 


0 tantj), „ . e 0 tancjK 

e o (mB - r^sin3 e o T ) 


From Fig. A. 2: 


Lp2 “ HE/2 = 1/2 r^cos3e^o tan< ^ 
From Equations A.24d and A.24e: 


P 0 L „ = 1/2 Y ricos 2 3e®o tan< ^(mB - x. sin3e^o tan< ^) 

L PZ S 1 

Determination of P -L , 

si si 

The average value of earth pressure on vertical plane EF is given by: 



!/2 K o . Ys *(D -D e ) 


(A. 23) 


(A. 24a) 


(A. 24b) 


(A. 24c) 


(A.24d) 


(A.24e) 


(A. 25) 


(A. 26a) 
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(A. 26b) 


The force P g ^ may be written as follows: 

p si - 1/2 VV® -V 2 

The moment arm L g ^ is given by: 

L sl ■» 1/3 (D + 2D b ) (A. 26c) 

Combining Equations A. 26b and A. 26c and making the substitutions D =* mB and 
= m'B yields the following expression: 

P sl L sl " 1/6 Y s ® 3 K o '<* " »■*)*<» + 2m') , (A. 27) 

Determination of P „L . 

sz s2 

The average value of earth pressure on the vertical plane OG is given 

by: 


p s2 - 1/2 

The force P g 2 may be written as: 


(A. 28a) 


P s2 - 1/2 R*Y s *D (A. 28b) 

The moment arm L g2 = D/3, consequently; 

P s2 L s2 “ 1/6 K *Y S *» 3 (A. 28c) 

or hf substituting D = mB: 

P s2 L s2 “ 1/6 K Y s (A. 29) 

It should be noted that the coefficients K q and K which appear in Equations A. 27 

and A. 29, respectively, are lateral earth pressure coefficients. In computing 

the .N values presented later in this appendix, they are assigned the value 

K = K = l-sin<j). 
o T 
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Solution of Moment Equation 

The moment equation (Equation A. 12b) may now be resolved by taking 

2 ’ ’ •- 

L b = j r, and by substituting the expressions given in Equations A. 18, 

A. 23, A. 25, A. 27, and A. 29: 

F b - K r i 2 

+ Y g r^ 2 sinf3 cos 2 8 e 

+ ^ Y g c°s^8 e ^o tai *==(mB-r^ sing e ^o tan< ^) 

+ J K Y s (B 3 /^) m 3 

- 7 K Y e (B 3 /r. ) (m-m 1 ) 2 (mf2m') (A.30) 

H O S 1 


Development of Expression for 

Because OC is a slip plane, stresses on plane OC can be represented by 
point A on Mohr's diagram. The plane OA is not a slip plane; however, 
stresses on this plane can be determined by the intersection of envelope 
T - atand and the Mohr's circle at point C. From the geometry of Mohr's 
circle, the following relationships can be written: 

OD = = 00 T +0'D (A. 31a) 

where 00 1 = (Q’A/sincJ>) and 0 'A = (AB/cos^)), and 



00’ * T~" . •: 

coscp sxn<p 

(A. 31b) 


O'D = O'C sin(2Y-4>) 

(A. 31c) 

where 

O'C = O' A = (AB/cos<t>) 


and 

q. d _ AB sin(2Y-4>) 

COS<f> 

(A.31d) 

Recognizing that AB = 

T^, and substituting Equations 

A. 31b and A.31d into 

Equation A. 31a: 



t r i 

^b 1 

O r . , + sin(2Y-<f>) 

p cos<}> sinifi 

(A. 32a) 
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By substituting = Octant}) into Equation A. 32a: 


a 

P 


V 


Normal forces on planes OC and OA 


F 

P 


a 

P 


[l+sintfr sin(2y-(j))] 
cos 2 <|> 

are given by: 

r and F, * a, r, 
o b b 1 


(A. 32b) 


(A. 33a) 


or 

V F b ■ VVW < A - 33b > 

Recalling the identity sina = cos (90° -a), the following relationship may 
be written: 


cos(y-<J>) 

r . r 

1 O costp 


(A. 33c) 


By substituting Equations A. 32b and A. 33c into Equation A. 33b 


F 

P 


, 1 + sinji sin(2y-4>) 

b cost}) cos(y-<|)) 


(A. 34) 


By substituting F^ from Equation A. 30 into Equation A. 34, the following 
relationship is obtained: 


F 

P 


1 + sintfr sin(2Y~<i>) 
cost}) cos(y-<}») 




, 1 2 jo Zo 30 tan<|) 

+ -j- Y g r -^ sinB cos 8 e o 

. 3 Zo 20 tantj) - _ . 0 0 tantj> s 

+ Y g cos Be o r (mB-r^sinBe o T ) 


+ | K Y g (B 3 /^) m 3 

" 4 K o Y s (B3/r i ) (nH-2m')| (A.35) 

Denoting the vertical component of F as F where F = F /cos6 and 

p pv pv p 

considering the vertical equilibrium of the force qB/2 exerted downward 
on the penetrometer and the resisting force, the following relationships 
are obtained: 
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(A. 36a) 


and 


qB = 2 F cos(^-6) ^ 


* -i*i**^ 


(A. 36b) 


Considering the weight (W) of the soil wedge for flat-ended penetrometers , 
Equation A. 36b is modified as follows : 


2 F 


£ cos OlJ-5) _ W 
B cos6 B 


(A. 36c) 


Also, the following relationship should be applied to the expression for F p : 


B cos (~Y— <f>) 


(A. 36d) 


12 cosiJj cos<j> 

Then substitution of Equation A.36d into Equation A. 35 yields the following: 

1 + sind) sin(2y-(|)) 


coscj) cos(y-(})) 

i Y B 2 I 

8 1 s cos 2 ^ cos 2 (f) 0 


. 1 _ 2 cos 2 (y-<|>) J0 2 o 36 tancj) 

+ ttY B z T 2a sin3 cos pe o 

16 's cos 2 i(j cos z cp 

3 n cos(Y-<f>) 2 d 20 tan<f >/ „ B 
+ 8 B co# cosf C ° S 6e ° r-2 


sinB e e o tan *) 
cosip coscp I 


+ i K J.costcosi m 3 

2 's cos(y-<p) 

1 b 2 C ° S f v C °f (»-’) 2 (nH-2m\) 

2 o 's cos(y-<p) 


(A.36e) 


Utilizing the relationship 


x- , cos(y-?<j>) j d 6_ta.u^ 

m' „ — v T/ - si n g e o , 

2 cosy cos<p 


Equation A.36e may be further simplified and combined: 
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„ 1 + sin<f> sin(2y-<i>) 

* 58 7 - 777 jl v “ 


COS COS (y-4>) 


— v •*tvS-* *^ss' _( Y~ ^ ) j 

8 cos 2 ^ cos 2 <J> 6 


+ 2 B 2 cos 2 e e 29 o ta "* (»-!„') 

8 ’s cosljj cos<j) v 3 

+ tKv b 2 s ° £L£g|| m 3 
2 s cos 


_ I K Y B 2 COSli* COS$ ( _ ml)2 (m+2int) 

2 o 's cos (Y-<t>) 


(A.36f ) 


Substituting the expression for F given in Equation A.36f into Equation 
A. 36c the following relationship is obtained: 


q = (Y g B) 


cos ($-<$) 


1 + sin<|> sin(2Y-4>) 

cos<S 


cos<j) cos (y-4>) 


1 cos 2 (-Y-(|>) T 
I 4 cos 2 ^ cos 2 <j) 0 

, 3 cos (y-<i)) in 20 tan<J> , 2 tx 

4* “7 — — : r J - ,, T l r cos 3 e o (m - m 1 ) 

4 cosy coscp N 3 


+ K 


COSip COS<P 3 

/ T\ ^ 

cos (Y"9) 


. K (m— m ’ ) 2 (m +2m') 

o cos (y-<j>) 


Y B tanij; 


(A. 37) 


Considering that q f = cN c + Y g B N , the following expression is obtained 


for the N factor: 

yq 
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cos 1 4- sintj> sin(2yH>) 

N yq * cos6 J cos<j> cos (Y~ 4*) 


[ 


1 cos ^ (Y-(j>) J 
4 cos*$ cos z 4> 0 


3 cos (y-6) 2 o 20 tan$ , 2 ,* 

+ -r ■ , ' ^ 4- cos 0 e o r (tn - -r m ) 

4 cosip cos<p -J 


... „ COSip COSfp __3 

+ K COS V« m 

- K C °°V 2 (B+ 2m'>l - ^ 

o cos Cy-4>) J 4 


where I fl is given by Equation A, 20. 


(A. 38) 
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CURVES FOR BEARING CAPACITY FACTORS 
FOR DIFFERENT VALUES OF ANGLE OF INTERNAL 
FRICTION (<{)) , BASE SEMIAPEX ANGLE (a) , 

BASE ROUGHNESS (6/<j)> , AND RELATIVE DEPTH OF BASE (D/B) 


Note: 


In computing values of N . the 

and K (see Equation A. 38) were 
o 


lateral earth pressure coefficients, 
assumed as K“K “l-sin<j>. 


K 
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FIG. A. 3 BEARING CAPACITY FACTOR N c FOR WEDGES (GENERAL SHEAF 
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ERMAL FRICTION, 4 > - DEGREES 


FACTOR N c FOR WEDGES (GENERAL SHEAR) 



























A COMPUTER PROGRAM FOR THE DETERMINATION OF BEARING 
BEARING CAPACITY FACTORS OF WEDGES AND CONES 


The following computer program which consists of a main program (NC or 

NGQ) and one subroutine (ANG) was written for the determination of bearing 

capacity factors N and N of wedges and cones, 
c Yfl 

Subroutine ANG calculates angle y (the topmost angle of the plane shear 

zone) from Equation 3.27. Programs NC and NGQ calculate bearing capacity 

factors N and N from Equations 4.5 and 4.9 respectively. Angle & is 
c yq 

calculated in the main program and Equations 4.26 and 4.32 were used for the 
shape factors in the calculation of cone bearing capacity factors. 


Data Input 

Control Card (I 10) 

Columns 1-10 Number of factors to be determined (N) 

Note: Each set of a, 6/<j> and D/B corresponds to a single bearing 

capacity factor to be determined. 


Penetrometer Property Card (3 F10.0) 

Columns 1-10 Semi-apex angle of penetrometer (a -ALFA) 

Columns 11-20 Roughness of penetrometer (6/<J>-FAS) 

Columns 21-30 Relative depth of penetrometer (D/B-DPT) 

N number of cards are required 

Notes: (1) Bearing capacity factors for both wedges and cones are 

calculated for 20°<<j><50® with increments of 5®. 

(2) Separate Input Data should be prepared for program NC 
and program NGQ. 
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LISTING OF THE COMPUTER PROGRAM FOR THE 
DETERMINATION OF BEARING CAPACITY 
FACTORS OF WEDGES AND CONES 
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KJ U U U U U U U 


PROGRAM Nr ( I NPUT ,OUTPi IT ) 


PROGRAM COME RES I STANCE Nf FACTORS 
N= NUMBER OF OAT A so ts 
A L F A s HALF APEX ANGLE 
FAS- ROUGHNESS 

UP T = DEPTH TO DIAMETER RATIO 
f i = angle of internal friction 

ATAR= SHAPE FACTOR 

DELT A= PEmFTromFT FR I O SOIL FRICTION AnG| f 
READ 1 , N 

1 FORMAT ( 110) 

C**-*-** READ AnD o-r In I' L'FDGE AmG|_E’ ooUGHN r S5» qFj . dFoTh *■**»■ **•#«.*• 
DO 30 A = 1 » N 
READ 1 0 »ALFA,FAS»DPT 
10 FORMAT CO Ft 0.0) 

PRINT pi , ALFA , FAS , DPT 

2i Format ( im >/,x »r-halfa=»fi •?, // 

1 5 X » 1 ohroughnf ss- ,f*-. ->,// 

2 AX * 1 0HREL. depths ,F« .-,//) 

ALFARsALFA/^7. ?fl 

F I = 2 0 • 

PRINT 24 

24 Format ( /•* i ox»?hf i •» i sx , - usps.i »i x ♦ *hg.a m a ♦ i /,x , .'h-iifTa, nx , b-hdflta 
1 HWLDGE FACTORS, 5X, IRHCONc FACTORS, /) 


DO ?5 J = 1 * 10 
F I R = F I / 5 7 . 2 O 

ATAR=i .0+0. ? + .( TAM (FIR) **a) 
q *** **£ ^ L.oUL T I ON oF ANr Lf ^A |V| /i'^^^***##'^**#'^*** 
CALL AN'u (FIs F a S » G A v A ) 

V6=SIN( FIR) /CuS(FIR) 

DELT A= FAS*FI 
DELTAR = DFLTA/S-7.2o 
IF (FA S - 1.0) 0 s , o ? , 3 "3 

3 2 F I S = 45 • — I T/2. 

IF ( ALFA-FIS) TS., t'S»1a 
16 P'S I = F I S+F I 


■K* >r a- -X- X- X- X- 


>r >v‘ > " -,r a' \ 


GO TO 17 
15 CONTINUE 
33 P S.I =90 • —ALFA 
17 P S I R = P S I / 5 7 , 2 8 

(2 -ic ~A? 'X* C ALGOL A i F ,,r T A «Y j T t~R a. T j ON ** **** * '• ” * *#***-***•# 

BE TA = 1 . *F I 
BE TAR-BETA/ 57. 2° 

0 A M A R = G AM A / 5 7 . 2 Q 
T E T A = 1 RO.-PSI +BETA-GAMA 
T E T A R = T E T A / 5 7 . 2 8 

uLUpi = o. 5 -;f ( Cqo { GAvAr-F IQ) ) * ( S.I n (biTAo) ) *lxp ( ( TETA**S I h (Flo) ) 
1 CoS (FIR) ) / ( CoS ( F I p ) -if COS ( p$ I P ) ) 

PRINT 44 » or DPT 

44 FORMAT ( / » 5 X , 4 HOC « = , F 1 n » 2 , / ) 

IF ( DFDPT-OPT) 11,11,12 
12 E T A = 1 8 0 » -P G I -G A "A 
SOTA=0. 

ET AR= ETA/5 7.23 

f 3=?. !fDP T » COo ( I I R ) * CQ~> ( PHI R ) / ( COS ( GA" 4R-F I R ) * ( EXP ( FTAR^Vr ) ) ) 
60R= ATAN ( T 1 ) 

BO= BOR* 5 7. 2 8 
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IF( BO-FI) 3 1 Q » 3 1 ft , 3 n 
311 BO=F I 
310 CONTI NUF 

DO 101 1 = 1 ,20 

T ? = ?. n 0 P T #' C 0 o i F I R ) «C0»<- PS I R 3 / ( COS ( GA'-AP-F I R } * { FX-o ( ( ( i o-O.-PSI 
1 -GAM A + BO ) / A 7 « 2 8 ) * V6 ) ) ) 

T3 = T?/( SQR T ( 1 . - T ? ** 2 ) ) ■ 

BN R- A TAN ( T" 5 ) 


B N = R N R 5 7.3 3 

I F ( A dS ( o,N- : 30 ) - d • i ) 3 0 ] » 3 0 1 * 3 0 ? 
r '01 GO TO 1.02 
302 dO= ( dN+BO ) / 2 • 

101 CONTINUE 

102 BFTA= BN 

dETAR- BETA/ 3 7 • PS 
PRINT 668,1 

668 F OR. 1 AT ( A X ? 1 OH I TER AT 1 0N= » I "> , / ) 

PRINT 6 6 ° , BETA 

6 6 Q F ORM A T ( /., X , 6 Hb<_ T A.* , F 7 • 1 » / ) 

C.**** CALcULA Tr WFnOF a No cON r p Ac TORS *#»*#*#■*■.* 
TLTA= l A0.-P61 HbETA-OAMA 
T E T A R = T E T A / r - 7 • ~ 3 
• GO TO 7 0 




11 CONTINUE 

L>1 = ( 7 ... XLQSl F IR) *CftS< PSIR )*UPT ) / < COS < GA v.Arj-F I o ) ■*C'"'G ( d f TAo ) *£yn [ TP" 
]R*V6) } 

OS I R = AT AN { D - 1. ) 

Q S I = QS I R 3-7 • 2 8 
o 0 T A = 0 • 

70 SOTaR=SOT A / 5 7 • 28 
C 1 = 3 I N ( F I R ) 

C?=COS( FIR) 

C3 = S I N ( 2 . *G AHAR-F I R ) 

C4 = EXP ( 2.*TFTAR*VA ) 

C-3 = S I N ( 7 • *SO TAR+F I R ) 

C6 = C0S ( 2 . -::-G AMAR-F J R ) 

C7=SIN(PSTR ) /C0S( PSIR) 

Cg = < l 1 • +Ci *C3 ) / ( C t 7 ) J < ( C^-C ) * ( i ... Tpff- >')/(( 0**9 ) ( Co-Va-: 

] (C3-Ci ) ) ) )*C 4 -i • ./Me. + ( Ca-s-Cp-s-La ) /O 
C9.= .{ V.6*C7*C6*< C3-Cl ) *Ca ) / ( C?* ( O’-Va * ( O -Cl ) ) ) 

ENC=C3+C9 
E N C O - A T A R «■ F N C 

PRINT ?a»FI , pSI , G A mA » tet a , delta, Enc , enco 

2 6 FORMAT i 4X ,f 1^ . 1 » /X , Flu . 1 , 3X , FlO . 1 , 9X , f lo , 1 , 3X , '“l 0 . 2 , 6..X ,rl A. 1 , 9X » c 
110.1) 

71 F I =F I +3 • 


25 CONTINUE 
30 CONTINUE 
END 
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SULiROUT I \)E A No ( F I , F AS > ON£W ) 

C CALCUUAT10M OF GAmA AnGl'F 

DFLT A = F AS*F I 
DELTAR=DELTA/57.?s 
IF ( FAS* tO. 1 . ) .GO TO ? 

I F ( FAS. EQ.O • ) GO TO 7 
Q S I = ( 1 • -F.AS ) * ( AS . +F I /? . ) 

QSIR-QSI A7.28 
F I R= F I / 57 • ? 8 

VALO= ( ol !\| ( uuL I AR ) /<_0S( ULlTAr ))*(]. +5 I N (FIR ) *5lN ( ? • *QS IR-F IR ) ) 
1— SI.M(FIR) *COS ( 2 • *OS I R-F I R ) 

DO 81 K = 1 » 150 
OS I =QS I +0 « 1 , 

QSIR=Q SI/87. 28 

VALN = ( S i R ( DFlTAr ) / CoS ( DElTAr ) ) » ( l . +SI m ( F I R ) *-S IN ( R • # 0 S I P— FIR) ) 
1 -SIN ( FIR ) *COS ( 2 • «-OS I R-F I R ) 

V AL = V ALN*VA LO 
1 F ( VAL ) 1 2 . >17 , 1 7 

11 VALO=VALN 
81 CONTINUF 
13 ONFW-OSI 
GO TO 4 

1 2 DOS I = ( VALM/ ( VAlm-V-ALO ) ) *0 . 1 
ON FW=QS I -DO S I 
GO TO 4 

3 QNEW=45 *+FI /?• 

GO TO 4 

2 QNEW = f). 

4 CONTINUF 
RETURN 
END 
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n n n n o n n 


PROGRAM NoO ( INPUT, 0!JTP-"T) 

c Cone, and weuul resistance mgo factor c 
N= NUMBER of BA T A A F T <?, 

A LF A = HALF APFX ANGLE 


FA As ROUGHNESS 

DPT = DEPTH TO DIAMETER RATIO 
FIs ANGLE OF I NTERNAL FR I FT I o« 

ALE M = SHAPE FACTOR 

D £ L T A = PENETROMETER TO SOIL FPJCTIPM A '“G L r 
READ 1 , M 

i format ( i io) 

c#**** read and print wedge angle » rough.vf c s» pfi , 


DO SO K = 1 , N 

READ ].0» ALFA, FAS, DPT 


10 FORMAT (SET 0.0) 


PRINT ? 1 > ALFA, FA F, DPT 

2 1 format t )H i , ax , <=.halfa=, fi o. ? , // 

1 GX » 1 ohroughnE fs= , Fg * ? , / / 

2 G X » i 0 HR E L • PE P T H = , r G • 7 , / / ) 


DFpT H ***#*%#■##'■#*##* 


ALFAR=ALFA/57 . 28 
PRINT 2 A 

24 format t / , i ux,2hfi , i sx ,ghps i , i x , ahgava * i ax a h t r t a , o x , g hd fl t a , r x , ? 

1 HWLDGL PACT ORG » GX, 1 2HCONE FACTORS, / ) 

F I = 2 0 , 


DO 2 5 J = 1 ,7 

F I R= F I / G 7 , 2 8 

A L LM = w • 6 + ( 1 « 5 / ( 1 • /D P I + ( 1 . r > / ( , *.+ ( T AN ( F I R ) -x- *- F ) ) ) ) ) 

C * * *• * o A L c U L A T I O N O F A N a L F G A ■'> A * * * * * * * * * * ' /! iv * * * * * * * * '* * * * 
CALL ANG ( F I , f AS , GAMA ) 

V6 - S I N ( FIR) /cost FIR ) 


■X- it if 


if -iS’ >7 if if if >f if if it 


DELTAS FAS*FI 
DELTAR=DELTA/g-7. 2 P 
IF (FAS-1.0) 77,72,83 

32 F-IS-A5.-FI/2. 

IF ( ALFA-FIS) IS, 15*1 A 

16 P S I = F I S + F I 
GO TO 17 

15 CONTINUE 

33 P 51=00 .-ALFA 

17 P S I R =P 5 I / G 7 . 2 B 

C*****CAl COLA r iOM BF A“tr L.F R c T A RY j T F R A T T ^'M*»***«-** »*******»*.« »**««*** 

beta=fi 

B E T A R = B E T A / 5 7 .28 
G AiMAR = G A '14/57.28 
TETA= LBO.-PSI+BETA-GaVvA 
T E T A.R = T E T A / 5 7 . 2 8 

DEDpT=to.s*COo< G Am A R -FIR >*SlN < PL TAR ) xLXP ( TETARx-Va ) ) / ( CO 5 < F I R ) *COS 
HP SIR)) ' 

PRINT 44 » B - B P T 

44 Format ( / , k a , /, hdc . - , f i > . 9 , / ) 

I F ( D.EDP T -DP T ) 3 1 , i i , ) 2 
12 DFDPT = DP T 

E TA= 180. -PS I -GAMA 
ETAR= FT A /G 7. 2 8 

} 1 = ( 2 .^DPTx-C0o ( F I R ) «COS ( PSIR ) ) / ( CD5 <GAhAP-Fi R ) *{ FXP ( FTAR*V*'> > > 

BOR = ATAi\ ITT) 

B 0 = BOR *5 7 . 28 
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IF (BO-FI) 310, 9)0, Mi 
311 30= FI 

310 CONTINUE 

DO 1.01 1 = 1,20 

T 2 = ( 3 . •;:-DPT-»(.Oi> ( Flo ) »C0S ( n$ Ip > ) / ( cos ( GA<'Ap-FI p ) ( LTD { ( < i <■> n*,_pc i _r 
1 A + BO ) / 5 7,28 ) *V6 ) ) ) 

T 8 = T2 / ( SORT ( 1 ,-T?**? ) ) 

BN R= A TAN ( T8 ) 

BN =B NR* 5 7. 88 

I F ( ABS ( BN-BU I-O.l) sO \ , 80 t , c, n ? 

R 01 60 TO 102 
802 BO= (BN+BO) /?. 

101 CONTINUE 

102 BETA= BN 

B E T A R = BETA/ 8 7,28 
PRINT 56 8, I 

868 •' ORMAT ( a X , i v-l i I ILQ AT I 0M.= » I 2 » / ) 

PRINT 669, BETA 
a 6 ° format ( a x , k h b l t a= , f 7 . i ,/) 

TETA= 1 8 O.-pSI-GAMA+utT a 
TFTAR= TET A/6 7. 28 
11 CONTINUE 

R AD I 5 -OF DP T / T A N ( 3 F T A R ) 

PRINT 48, RAD I 5 

4 5 FORMAT ( / , 8X »4HRD. - ,Fl u. 2 > / ) 

C *** **-CALCvJI_aT I ON of NOU F ACTutS £>•**•*’*■# ##«•*•» * * * * * 8-:,-x **•« **-■» *** 

W E I G H = S I N ( P 6 I N ) / ( 4 . * COM PSI R) ) 

I r ( AL r A-o 0 . ) 13 , 14,18 

13 WE I GH=0 • 

14 CONTINUE 

V1 = FXP ( 3 . *T F TARM/6 1 
V 2 = V 1 *CO$ ( BETAR )-C03 ( TFT AR-UE TAR ) 

V8 = V l *3 1 N ( BE TAR ) +8 I N ( TlTAR-.h_TAR ) 

V4=3.*V6^V2+V3 
V 5 = V 4 / ( 1 » +0 ,■'!■( V A - "' 2 ) ) 

U 1 = C 0 S ( 6 AM A R - r I R ) 

U2=COS( FIR) 

U 3 -CO 81 P S T R ) 

U4=COS ( c.FjAR ) 

u 8 = ( *i , / 4 . ) -x ( U 1 "• •* 9 ) :s V 8 / ( ( i 1 8 X 8 ) -;r- ( j 8 * ) j 

U7 = EXP ( 2. ::-TFTAR;:-Va ) 

ua=ui* (U4**? I/O n*u ? ) 

U9=(3./4. >*Ufl#U7*(DPTM 2./3. )*-^FDPT) 

D6=l .-SIN (FIR) 

02= ( ( D6*U8*u2 ) / ( Ji ) ) «•( ( DpT-DFDdT )*#;>) * ( Dr> 7 +•■>.* D c v "T ) 

D4 = ( ( D6*U9«-'J2 ) * ( DP T* *3 ) ) /!) 1 

Pi = ( U 3 + U q— D 2 + 1 ' 4 )*( C0$ ( PS I R-DE|_TA D ) /C^S ( NE[ TAo ) ) 

P2= ( ( i «+3I M ( F I R ) -:«• 81 N ( 2 « *GA|VAo-F Ip) ) / ( U?*.Ui ) ) •“ °1 
t M 6 A Q = P 2 - w E I G H 
E N G A C = A\ L E M *' £ N G A Q 

PRINT 2 6 5 F i ,p5l , G AmA , TL T M ,i>LL ^ A , L'.jSiAO , L 'iGA C 
2 1> FORMAT FAX, Fl^.l, /X»Fl'8,l,3X»Fi . ] , y X , F i j . 1 , 3 X • F 1 0 » 2 , 6 , F 1 , j. • 9 X • 
1F1 0. 1 ) 

FI=FI+5.0 
25 CONTINUE 
30 CONTINUE 
END 
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S U a R 0 U T I N £ A N G ( F I , F A S » D .« F w ) 

C calculation OF GA-'A ANGLE 

DELTA = F AS#F I 

L> L L T A R = D L L T A / 5 7 • 7 <3 

I F ( FAF.LQ.1 « ) GO TO ? 

IF (FAS. CO. 0. ) 'GO TO 
051 = ( 1 . -F AS ! * ( AS . +-F I / ? . ) 

GSIR-05I/57.28 
F I R = FI / 5 7 . 2 8 

valos? < si m ( dll Tar y/c n.*> ( ullTAr y ) * ( i . +s lm ( F i o y #s tr ( ? » #o5 i r-f ip)) 
1 -5 I N I F I R ) *.0 OS ( 7 . *o s I R-f I - ) 

DO 81 iv=1 , ) 50 
wS I =QS I +0 • 1 . 

OS l R=OS 1/57.78 

VA|_N = ( o I a.i (Dll i HP ) /v.OS ( w>L|_ T nr> ) ! -v ( i . +S I m ( F I n 5 **> I m ( o • *o$I r>_F I P > .) 
1-5 IN (FI E ) * CO S ( ? . * n 8 I R - f I li ) 

V AL = V A L N * V A LO 
IF (VAL) 12,18,11 
11 VALO-VALN 
81 COOT 1 NO F 
13 QNF''!=Q S I 
GO TO A 

1 7 Dro I - ( V Al N / ( VmL'I— VA i n) ) i 
UNFU=QS I -DOS I 
GO TO A 

3 ON£W=45.+FI /?• 

GO TO* 4 

2 QMFW=0. 

4 CONTINUE 
Re- TURN 
END 
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APPENDIX B 


SPREADER BOX CALIBRATION 

Pluvial compaction, the deposition of air-dry sand by allowing it to fall 
as a rain to build up a bed, was used to produce a large range of densities. 

Sand was stored in 55 gallon (~0.22m 3 ) drums between tests. Prior to deposi- 
tion, the frame travel micro-switches were adjusted and the spreader box was 
filled with sand using a barrel turner attached to a fork lift truck. Motion 
of the frame was iniated one foot from where the roller started to eject sand 
to ensure that the frame reached its required speed and that transient vib- 
rations had died out. A second micro-switch started the deposition which 
continued until a third switch was contacted stopping the rain of sand. A 
fourth micro-switch reversed the direction of travel allowing the frame to 
move back and forth spreading out about a 1.0 in. (25.4 mm) thick sand layer. 

The spreader box was intermittently elevated to maintain a constant height drop. 
The procedure was repeated until the required sample height was attained. 

In order to determine the uniformity of the sand layer deposited by the 
spreader box, the sand density at various locations in the layer was measured 
by weighing the quantity of sand deposited into volume calibrated cylinders. 
Cylinders of 3.0 in. (76.2 mm) diameter and 3.0 in. (76.2 mm) height were 
used as recommended by Kolbuzewski (1948). In several tests, cylinders were 
placed along the transverse and longitudinal direction of the box motion. 

Figs. B. 1 and B.2 show the variation of sand density along these directions 
of the box. In several tests, cylinders were placed at different elevations 
to check if there was any density change due to additional sand deposition. 

It was found that the variation in density was negligible, even for loose sand 
deposits. The overall average density was determined by weighing the quantity 
of sand deposited into the test box. 

The density of sand deposited using the spreader box is related to 
roller speed for a given gap width and height of fall so that once cali- 
brated, it is possible to reproduce any required density by properly adjusting 
the roller speed. The calibration curve obtained for Monterey Sand No. 0 
is shown in Fig. 8.3. It may be seen from Fig. B.3 that obtainable relative 
density range is quite large, 34 to 100 percent. 
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